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Abstract: In order to calculate QED corrections to hadronic physical quantities by means 
of lattice simulations, a coherent description of electrically-charged states in finite volume 
is needed. In the usual periodic setup. Gauss’s law and large gauge transformations forbid 
the propagation of electrically-charged states. A possible solution to this problem, which 
does not violate the axioms of local quantum field theory, has been proposed by Wiese and 
Polley, and is based on the use of C* boundary conditions. We present a thorough analysis 
of the properties and symmetries of QED in isolation and QED coupled to QCD, with 
C* boundary conditions. In particular we learn that a certain class of electrically-charged 
states can be constructed in this setup in a fully consistent fashion, without relying on gauge 
fixing. We argue that this class of states covers most of the interesting phenomenological 
applications in the framework of numerical simulations. We also calculate finite-volume 
corrections to the mass of stable charged particles and show that these are much smaller 
than in non-local formulations of QED. 
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1 Introduction 


Electromagnetic interactions contribute at the order of a few percentage points to masses, 
decay rates and scattering cross-sections of hadrons. Nevertheless these small contributions 
cannot be ignored if one is interested in quantifying isospin breaking effects like the charged- 
neutral mass splittings of baryons and mesons, or when one aims at percent accuracy in the 
calculation of hadronic matrix elements. In these cases first-principle theoretical predictions 
can be obtained only by means of lattice techniques, which require a consistent formulation 
of QCD-I-QED in finite volume. 
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The problem addressed in this paper arises every time one needs to produce an electrically- 
charged state in a hnite periodic box, as for instance in the calculation of the proton mass, 
and is intrinsically related to the dynamics of the zero-modes of the gauge field. In a torus 
with periodic boundary conditions for the gauge fields. Gauss’s law implies that only neutral 
states belong to the physical Hilbert space of the theory. One might think to overcome this 
limitation by gauge-fixing. For instance in Coulomb gauge the Gauss’s law is locally solved 
and the Hilbert space splits in sectors labeled by the total electric charge. However states 
generated by electrically-charged local operators in Coulomb gauge are also charged under 
large gauge transformations which survive a local gauge-fixing procedure. Because of this, 
even after gauge-hxing, the two-point function (^/)(x)'0(y)) vanishes if x and y are separated 
in a periodic box. In practice large gauge transformations act on the gauge held by shifting 
the global zero-modes Therefore the obstructions to the propagation of 

charged particles on a periodic torus can be traced back to the functional integration over 
the global zero-modes. 

A possible solution to this problem can be found in ref. [1] where the hrst lattice calculation 
of the electromagnetic mass splitting of nucleons and light pseudoscalar mesons has been 
attempted. The proposed solution consists in quenching a particular set of Fourier modes of 
the gauge held, in such a way that the global zero-modes decouple from the dynamics. A lot 
of theoretical and algorithmic progress has been made after the pioneering work of ref. [1], 
particularly in the past few years, leading to recent determinations of the electromagnetic 
mass splitting of light pseudoscalar mesons and light baryons, see refs. [2-9] for recent 
works on the subject (see also ref. [10] for the discussion of a method to calculate QED 
radiative corrections to the leptonic decays of pseudoscalar mesons). All these works rely 
on hnite-volume formulations of QED obtained by quenching some Fourier modes of the 
gauge held.^ 

The particular formulation called QEDl is obtained by quenching the spatial zero-modes of 
the gauge held at any time, i.e. by enforcing the constraint A^{t, 0 ) = d^x x) = 0. 
As opposed to other formulations, QEDl has a well dehned transfer matrix. However the 
constraint 0) = 0 is non-local. Even though one can argue that the modihcation 

generated by the constraint is a hnite-volume effect, many properties of local quantum held 
theories are not automatically guaranteed for QEDl- Among these we mention renormal- 
izability, volume-independence of renormalization constants, the validity of the operator 
product expansion and of the Symanzik improvement program. Mild violations of locality 
may preserve some of these properties but this needs to be shown explicitly case by case. 

QEDl has been studied at one-loop in perturbation theory in refs. [2, 4, 14]. The quench¬ 
ing of the zero-modes does not generate ultraviolet divergences at one loop, other than 
the inhnite-volume ones. However it does generate unusual phenomena, for instance par- 

^ Recently other approaches have been proposed. In ref. [11] the zero modes of the gauge field are lifted 
by adding a mass term for the photon. The proposal of refs. [12, 13[ consists in combining QCD matrix 
elements extracted from finite volume simulations with infinite volume QED kernels. 
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tides and antipartides do not decouple in the non-relativistic limit [15, 16]. This can be 
seen as a failure of the effective-theory description which is not surprising if the underlying 
microscopic theory is non-local. On the other hand the numerical results of lattice simu¬ 
lations of QEDl performed in refs. [2, 5] might be viewed as reassuring evidence that the 
non-localities of QEDl have only mild effects on the hadronic spectrum. Nevertheless we 
believe that QEDl is not sufficiently well understood at all orders in perturbation theory. 
Our approach is to eliminate any potential problems at the root, by seeking a consistent 
formulation of the finite volume theory that does not require quenching dynamical degrees 
of freedom. 

In this paper we consider a local solution to the problem of charged particles in finite 
volume. This solution is not new, it has been proposed in [17-20] and consists in enforcing 
C* boundary conditions for all fields along the spatial directions, i.e. in requiring that the 
fields are periodic up to charge conjugation. In this theory, which we refer to as QEDc, the 
zero-modes of the gauge field are absent by construction because A^{x) is anti-periodic in 
space, and the classical problems of the periodic setup are avoided from the very beginning. 
We show that a complete description of a certain class of electrically-charged states can 
be obtained without relying either on perturbation theory or on gauge-fixing. As we shall 
discuss in detail, this class of states covers most of the relevant spectroscopic applications 
and includes the proton, the neutron, the charged pions, the charged kaons, the charged D 
and B mesons and the baryons. The proposed construction is based on the fact that 
C* boundary conditions break the global gauge symmetry group U(l) down to its discrete 
subgroup Z 2 . In other words charge conservation is partially violated by the boundary 
conditions. The full group of gauge transformations splits in two disconnected components: 
the subgroup of local gauge transformations which are connected to the identity, and the set 
resulting by the composition of local gauge transformations with the nontrivial global gauge 
transformation. In this setup one can construct states that are invariant under local gauge 
transformations but not under global gauge transformations, and these can be identified as 
electrically-charged states. 

Along with charge conservation, C* boundary conditions partially violate flavour conserva¬ 
tion. This happens because flavour-charged particles traveling once around the torus turn 
into their antiparticles, and therefore change their flavour content. Being associated with 
the propagation of massive colorless particles, these effects are exponentially suppressed 
with the volume. We study in detail the pattern of flavour violation in QEDc, particu¬ 
larly in the case when electromagnetic interactions are coupled to QCD, and quantify these 
effects in the framework of a generic effective theory of hadrons. In particular we show 
that, although the and baryons can mix with lighter states because of the boundary 
conditions, the exponential suppression is so strong that these mixings can hardly represent 
a problem in numerical simulations. 

Finite-volume effects on the masses of charged particles are considerably smaller in QEDc 
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than in QEDl- When these corrections are expanded in a power series in 1/-L, at order 
aem in both theories the 1/L and 1/L^ finite-volume corrections to the mass of a charged 
particle are universal, i.e. they do not depend on the spin and on the internal structure of 
the particle (for QEDl see refs. [2, 4, 15, 16, 21]). We show that these universal corrections 
are always appreciably larger in QEDl than in QEDc. For instance at rriL = 4 we gain a 
factor of about 2 with three C*-periodic spatial directions and a factor of about 5 with a 
single C*-periodic spatial direction, see figure 4. More importantly, the spin and structure- 
dependent corrections are 0{1/L^) in QEDl, while they are only 0{1/L^) in QEDq. This 
extra suppression can be seen as a direct effect of locality. 

The paper is organized as follows. In section 2 we introduce C* boundary conditions and 
study the symmetries of QEDc. In section 3 we introduce the gauge invariant interpolating 
operators for charged particles and study their properties. In section 4 we couple electro¬ 
magnetic and strong interactions and study the symmetries of QCD-|-QEDc. In section 5 
we discuss the finite volume corrections to the masses of charged hadrons. In section 6 
we discuss the details of the lattice implementation of C* boundary conditions and of the 
proposed gauge invariant interpolating operators. We draw our conclusions in section 7. 

The paper contains four appendices with the explicit derivation of some of the results 
presented in the main body of the paper. The material discussed in the appendices is 
technical and some of it is, we believe, original. Appendix A presents a detailed study 
of some flavour-violation processes in QCD-|-QEDc, in the context of a generic effective 
theory of hadrons. This analysis requires an extension of the techniques developed to study 
finite-volume effects in [22], and it is complicated by the need to keep track of flavour flow 
and violations through all possible Feynman diagrams. In appendix B we give an ab-initio 
derivation (i.e. without using an effective description of hadrons) of the power-law finite- 
volume corrections on the mass of charged hadrons in QCD-|-QEDc. The coefficients of the 
expansion in powers of 1/L are expressed in terms of physical quantities, i.e. derivatives of 
the forward Compton amplitude for the scattering of a soft photon on the charged hadron. 
The authors are convinced that the technology developed in these appendices will find other 
uses in the field. 


2 QEDc 


In this section we introduce the finite-volume theory QEDc and study its symmetries. 
For simplicity, we consider the case of a maximally symmetric torus with linear size equal 
to L, with fields obeying C* boundary conditions in all space directions. The Euclidean 
time direction can be either infinite or compact with linear size T. In the latter case the 
corresponding boundary conditions for the fields will be left unspecified. Common choices 
are periodic, Schrodinger Functional (SF), open or open-SF boundary conditions. 
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The action of QEDc is given by 


S[AM 



The held strength and covariant derivative are dehned as 


( 2 . 1 ) 


F^,u{x) = d^A^{x) - dyA^{x) , 

Dl = d^-iqfA^, ( 2 . 2 ) 

^ .. —^ ■<— 

where the left-right derivative — dy) is dehned in terms of the partial derivative 

acting to the right and the partial derivative acting to the left. In our notation qj 
is the electric charge of the /-th havour normalized to the electric charge of the positron 
(i.e. qf does not include the coupling constant e). Throughout the paper we use this 
normalization for the electric charge. 


Fields obey C* boundary conditions under translations in the three space directions, 


A^{x + Li) = ^^(x) = -Af,{x) , 
ipf{x + Li) = V’/(x) = C"(x), 

tjjfix + Li) = iijix) =-'ipJ{x)C , (2.3) 


where Li is L times the unit vector in direction i. The charge conjugation matrix can be 
taken to be any invertible matrix C with unit determinant such that 

= - 7 j , (2.4) 

where 7 ^ are the Euclidean gamma matrices. In four dimensions such a matrix exists and 
satishes 

= -C , = C-^ , (2.5) 

independently of the particular representation of the gamma matrices. 

Notice that the action density eq. (2.1) is the same as in inhnite volume and it is therefore 
invariant under charge conjugation. Since a shift of a period in space corresponds to charge 
conjugation, the action density is periodic in space. 

We are now going to study the symmetries of QEDc, in turn gauge transformations, spatial 
translations, parity and flavour symmetries. 
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2.1 Gauge transformations 


Gauge transformations are defined in the usual way 
= A^{x) + d^a{x) , 

= . (2.6) 

Only gauge transformations that do not change the boundary conditions of the fields are 
admissible. Translating the transformed field by a period along a spatial direction yields 

+ ^i) ~ + Li) + dfj_a{x + Li) 

= -An{x) + d^a{x + Li) = + dfj_[a{x + Li) + a(x)] . (2.7) 

The transformed field ^ (x) is anti-periodic if and only if the gauge transformation satisfies 

dfj,a{x + Li) = -d^a{x) , ( 2 . 8 ) 

i.e. a(x) can be decomposed into an anti-periodic function plus a generic constant. The 
boundary conditions for fermions constrain this constant. Translating a fermion field by a 
period along a spatial direction yields 

V;[“l(x + Li) = + Li) 

= . (2.9) 

The transformed field ^|;^j)‘\x) satisfies C* boundary conditions if and only if an integer Uf 
exists such that 

a{x) = I3{x)+ , I3{x + Li) =-I3{x) . (2.10) 

Notice that this equation has to be satisfied for all fermion fields and for any pair of charges. 
In the physically relevant case^ all charges Qf are integer multiples of an elementary charge 
Qei, therefore the gauge transformation a(x) preserves the boundary conditions of all fields 
if and only if an integer n exists such that 

a{x) = I3{x)+ — . (2T1) 

Qel 

Quantization of the electric charge can be seen as a consequence of the fact that the gauge 
group is the compact U(l). A generic gauge transformation is assigned by choosing a 
phase factor A(x) = in each point of spacetime. A matter field with charge qf 

transforms with A{x)'^f where qj = qf/qd is an integer, i.e. accordingly to some irreducible 

two of the charges have irrational ratio, then one of the n/ has to be zero and consequently a{x) has 
to be anti-periodic 
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representation of the gauge group U(l). This analysis can be restated in terms of operators: 
given the electric-charge operator Q, the generator of global gauge transformations is 

Q = —, (2.12) 

Qel 

and has only integer eigenvalues. C* boundary conditions break the U(l) group of global 
gauge transformations. In fact eq. (2.11) implies that the only allowed global gauge trans¬ 
formations are A = ±1, i.e. the global U(l) is broken down to Z 2 . Breaking of the global 
U(l) implies a partial violation in electric-charge conservation: Q is not conserved but the 
quantum number ( — 1)*^ is. The origin and consequences of this violation will be discusses 
in more details in subsection 2.4 for the case of QEDc in isolation, and in section 4 for the 
case of QCD-|-QEDc. 

Eq. (2.11) implies that the group of gauge transformations is disconnected. Only gauge 
transformations with n = 0 , i.e. with a{x) anti-periodic in space, are continuously con¬ 
nected to the identity. We will refer to these gauge transformations as local gauge transfor¬ 
mations. Note that the large gauge transformations have a very simple structure (they are 
just the composition of a global gauge transformation and a local gauge transformation). 
This contrasts with the case of periodic boundary conditions in space, where large gauge 
transformations are linear in the coordinates (i.e. a{x) = ‘I'KnXijL with some integer n). 


2.2 Translations 


C* boundary conditions preserve translational invariance and charge conjugation. Even 
though in inhnite volume the momentum and the C quantum number are unrelated, this 
is not true in QEDq. Eqs. (2.3) imply that the translation of a generic (elementary or 
composite) held (t){x) by Lj is equivalent to a charge conjugation 

(l){x + Li) = (jf {x) . (2.13) 


The C-even and C-odd components of the held 4){x) are 


0±(x) 


(f{x) ± (ffix) 


(2.14) 


(/>+(x) is periodic in space while 4>-{x) is anti-periodic. The two components have different 
Fourier representations. Since we want to leave the time boundary conditions unspecihed, 
we expand our helds in the time-momentum representation, 

(f>±{x) = ^ '^±(^o,p)e*P'' , 

p6n± 
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(2.15) 


where 11+ is the set of periodic momenta and n_ is the set of anti-periodic momenta, 


n+ = l^n I n e z^l , 

n_ = 1^ (2n -I- h) I n G Z^, h = (1,1,1)| . (2.16) 

Notice that the Afj_{x) held is C-odd and it has only the anti-periodic component, while the 
helds '4’fix) contain both, 

M^) = p E E(^o,p)e*P" , 

pen_ 

V’/,±(a;) = E E±(®o,p)e*P^ . (2.17) 

p6n± 

The two V'/,± components of the fermion helds satisfy the (anti) Majorana condition 

V’/,±(a;) = ±C-^['4f,±f (x) . (2.18) 

2.3 Parity 


Even though not in a trivial fashion, parity is conserved by C* boundary conditions. Under 
parity the helds transform like 

Ao(x) Ao{xp) , 'ipfix) ??p7oV'/(a^p) , 

Akix) -Akixp) , 4^fix) r]*p4!f{xp)'yo , (2.19) 

where xp = (xq, — x). In inhnite volume rjp is a generic complex phase. For each choice of 
rjp one obtains a different but equally good parity operator. A customary choice amounts 
to r]p = 1. However the parity operator dehned in this way does not commute with the 
charge conjugation operator that we have used to dehne the C* boundary conditions. A 
more natural choice is rjp = i. The corresponding parity transformation V commutes with 
the charge conjugation operator. This can be shown explicitly by acting on the elementary 
helds with charge conjugation C hrst and parity V after, and by comparing the result with 
the same operations applied in reversed order. For example, in the case of the fermion held 
we have 

V’/(x) C~^'ipjix) A -iC~^jQipJixp) , (2.20) 

and 


'ipfix) A z 7 oV'/(xp) A z 7 oC' ^'ipjixp) . 


( 2 . 21 ) 



The results of the two transformations are shown to be equal by using C ^ 7 (fC' = — 70 - 
The reader can check that this conclusions applies to the other fields. 

Since V leaves the action and the C* boundary conditions unchanged, it is an exact sym¬ 
metry in finite volume. Even though parity will play no special role in this paper, we notice 
that the parity transformations can be easily used to construct operators that have dehnite 
parity. 


2.4 Flavour symmetries 


C* boundary conditions violate flavour (and consequently electric charge) conservation. 
The violation arises because a flavour-charged particle flips the sign of its flavour content 
by turning into its antiparticle when it travels once around the torus. We are now going to 
show that flavour is violated by two units at the time in this process and that this effect is 
exponentially suppressed with the volume. In this subsection and in section 4 we will argue 
that flavour violation does not represent a limitation to the use of C* boundary conditions 
in most of the relevant applications. 

We start by considering the theory with a single species of charged particles with unit 
charge, e.g. the electron. In this case flavour coincides with the electric charge Q and with 
the generator Q of global gauge transformations. The detailed way charge conservation is 
violated by finite-volume effects can be easily understood by means of Feynman diagrams. 
We assume here some gauge fixing that we do not need to specify at this level. The theory in 
hnite volume has the same interaction vertex as the infinite-volume one which, in particular, 
conserves electric charge. The violation of charge conservation is visible in those terms in 
the action that are sensitive to the C* boundary conditions, i.e. the ones containing spatial 
derivatives. In other words, charge violation is generated by the propagators, which we will 
discuss in detail. 

In order to write down the free propagators, one needs to keep into account the fact that 
a free particle is able to travel around the torus. If it travels once around a direction with 
C* boundary conditions, the particle turns into its antiparticle. The winding numbers of 
the particle world-line around each spatial direction can be organized into a vector n e Z3. 
By dehning 

3 

(n) = ^^rii mod 2 ( 2 . 22 ) 

i=l 

we can separate those winding numbers characterised by (n) = 1 that flip the electric charge 
of the particle from the winding numbers characterised by (n) = 0 that do not. We do not 
need the explicit expression of the gauge field propagator as the photon carries neither 
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electric nor flavour charge. Concerning the matter field, in coordinate space we have 


{'ip{x)ip{y)) = 

X 2/ _ 

->- - 

S{x -y + Lim) , 

(n)=0 

(2.23) 

{i;{x)'ip'^ (y)) = 

II 

- S{x - y + Lini)C~^ , 

(n> = l 

(2.24) 

{x)'ip{y)) = 

II 

yy CS{x -y + Lim) , 

(n> = l 

(2.25) 


where S{x) is the infinite-volume fermion propagator. Notice that the and 

propagators vanish in infinite spatial volume as the sums in eqs. (2.24) and (2.25) do not 
include n = 0. They are precisely the source of violation of charge conservation. The 
violation is not arbitrary, but amounts to a AQ = ±2 every time one of these propagator is 
inserted. This shows explicitly that the electric charge Q is not conserved, but the quantum 
number ( — 1)*^ is, which means 

AQ = 0 mod 2 . (2.26) 

Time evolution mixes all sectors with odd electric charge among each other, and all sectors 
with even electric charge among each other. For example a single-electron state can mix 
with a three-electron state but not with the vacuum, see figure 1. This in particular means 
that, chosen some suitable interpolating operator as we will discuss in section 3, single¬ 
electron states can be selected by looking at the leading decaying exponential in two-point 
functions. However two-electron states cannot be extracted in the same way, as the leading 
decaying exponential in a two-point function constructed with an operator with charge 
equal to 2 will select always the vacuum. As we will discuss in section 4 this is sufficient in 
most of the interesting low-energy applications in QCD-|-QEDc. 

From the discussion above it is obvious that the violation arises only from the charged 
particle that travels at least once around the torus. If the fermion is massive we have 

~ {x)i){y)) ~ " 6“”"^ , (2.27) 

for L —)• oo, and charge-violating diagrams are exponentially suppressed. 

In the case of Aj flavours the infinite-volume theory has a U(l)'^./' flavour symmetry cor¬ 
responding to independent phase rotations of each flavour. We will denote the generator 
of the /-th U(l) by Ff. Notice that the electric charge is a linear combination of the 
flavour-symmetry generators, 

Nf 

Q = ^qfFf, (2.28) 

/=i 
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Figure 1: (a) Diagram contributing to the e+e'*' —)• 7 process, which involves one e"*" 

traveling around the torus and flipping charge, (b) Diagram contributing to the —>■ 

e“ process. 

where (//is the electric charge of the /-th flavour. In infinite volume each Ff is conserved 
independently. C* boundary conditions break the flavour symmetry group down to a 
and this implies that only each ( — 1 )^/ is conserved, i.e. violations can occur only in 
multiples of two, 

AFf = 0 mod 2 . (2.29) 

Notice that the AQ has to be a multiple of 2 only if all flavours have equal electric charge. 
In the general case it is replaced by eqs. (2.26), (2.28) and (2.29). This observation will 
play an important role in section 4 where we will discuss QCD coupled to QED. 


3 Gauge-invariant interpolating operators 


We are concerned with physical observables, i.e. observables that are invariant under local 
gauge transformations. Often these observables are extracted from intermediate quantities 
defined in a particular gauge. For instance masses of charged particles are usually extracted 
from the long-distance behaviour of two-point functions after the photon field has been 
gauge-fixed. Although this is a necessary step in perturbation theory it can be completely 
avoided non-perturbatively without adding any particular complication. Keeping in mind 
that no issue arises with gauge fixing for QED and that a particular gauge can be chosen 
at any time, we think that it is more natural to rely on a completely gauge-invariant 
formulation. In this section we show how to construct states that are invariant under local 
gauge transformations and electrically charged at the same time, i.e. they have (— 1 )*^ = — 1 . 

® If n/ out of the Nf flavours are degenerate (i.e. same mass and same electric charge), the U(l)"t 
flavour subgroup is lifted to a U(n/) flavour symmetry. C* boundary conditions break this down to its 
natural 0(n/) subgroup. We mention this special case for completeness, but it is not relevant for the 
purpose of this paper. 
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This will be achieved by acting with suitably-constructed interpolating operators on the 
vacuum. Even though we discuss primarily how to apply this construction to the calculation 
of charged-particle masses from two-point functions, the same interpolating operators can be 
used to extract other physical quantities, e.g. decay rates, in a completely gauge-invariant 
fashion. 


To simplify the notation in this section we consider a single matter held with charge q. 
The generalization of the following discussion to the case of several havours with different 
charges is completely straightforward. Consider the operator [23] 

where is the matter held and J^{x) is a generic function or distribution that satishes 

di,Jfj,{x) = -F U) = -J^,{x) . (3.2) 

In case of periodic boundary conditions in time J^{x) is chosen to be periodic as well. Under 
a global transformation '!/)(x) —)■ the above operator transforms like 'I'j(x) —)• 

e*'?"Tj(x), which implies that in inhnite volume would have electric charge equal to 

q. In hnite volume we have already noticed that a can be only 0 or vr/g, which implies that 
the operator 'I'j(x) has quantum number (—1)*^ = —1. The non-local factor 

©(x) = A^{y)j^{y-x) (3 3^ 


transforms under a local gauge transformation that is anti-periodic in space as 
0(x) —^0(x) d.S oi(y)df,Jf,{y-x) 

= 0(x) . (3.4) 

Notice that the product a{y)J^{y—x) is periodic with respect to y. Given also the boundary 
conditions in time, no boundary terms arise from the integration by parts. The extra factor 
^-iqa{x) obtained by gauge-transforming 0(x) cancels the analogous factor coming from the 
transformation of “ipix), making 'kj(x) invariant. 


Summarising, the non-local operator Tj(x) has (—1)^ = — 1 and is invariant under local 
gauge transformations. It also satisfies the same boundary conditions as the held ip{x)^ and 
therefore operators with dehnite momentum can be easily constructed by considering the 
C-even and C-odd components of 'kj(x) as done in eqs. (2.14) for a generic operator <j){x). 


If the function J^{x) is chosen to be proportional to (5(xo), then the operator 'I'j(x) is 
local in time, i.e. it is a function of the elementary helds at the time xq only. In this 
case 'I'j(x) maps naturally to an operator acting on the Hilbert space. The state 'l'j(x)|0) 
obtained acting with the interpolating operator on the vacuum is invariant under local gauge 
transformations and has electric charge (—1)*^ = —I- By decomposing the Euclidean two- 
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point function ('I'j(x)^j(O)) in decaying exponentials in xq, one can extract the spectrum 
of the gauge-invariant Hamiltonian. The energy levels are gauge-invariant by construction 
and they do not depend on the particular choice of J^{x), as they are a property of the 
Hamiltonian rather than of the interpolating operator (as long as this is local in time). 
We will refer to the energy of the lightest state propagating in the Euclidean two-point 
function as the finite-volume mass of the charged particle. We assume that this quantity 
has an infinite-volume limit which can be interpreted as the mass of the charged particle.^ 

The whole construction presented above is based on the assumption that solutions of 
eq. (3.2) exist. If periodic boundary conditions were employed in all spatial directions 
eq. (3.2) would have no solutions. In the case of C* boundary conditions we will construct 
explicitly some possible choices for the function J^{x). The first one is defined by the 
equations 

Jo(x) = 0 , Jfc(x) = (I(xo)5fc^>(x) , dkdk^{-K) = (5^(x) , (3.5) 

where x = (xo,x) and <h(x) is anti-periodic. An explicit (convergent) representation for 
<h(x) is given in terms of the heat-kernel 

1 poo 

$(x) = -js du g-«p2+mx (3 6) 

pen_ 

With this choice the operator 'I'j(x) can be written like 

^^(x) = ^(^) ^ (3 7) 

Notice that in Coulomb gauge 'I'c(3^) = '>p{x), and therefore the gauge invariant correlator 
(\I'c(x)Tc(y)) is identical to usual correlator {'ip(x)'ijj{y)) in Coulomb gauge. In other words, 
'kc(a^) is the unique gauge-invariant extension of the operator 'ip(x) defined in Coulomb 
gauge. This in particular shows explicitly the gauge-invariance of the mass extracted in 
Coulomb gauge. 


Another possible choice is given by 
sgn(xfc) 5(a 


(3.8) 




Once this equation is inserted in eq. (3.1), it yields the following interpolating operator 

Tg(a:) = , (3.9) 


This choice generates a string wrapping around the torus along the direction fe, chosen 
among the ones with C* boundary conditions (see figure 2). The operator 'I's(x) is less 

^This issue is not trivial in QED because of the absence of a mass gap. See for instance the discussion 
in chapter 6 of [24] or chapter 6 of [25], and references therein. 
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Figure 2: Graphical representation of the interpolating operator 'I's defined in eq. (3.9). 
The black circle represents the electric charge, and the white circles are the image anti¬ 
charges. The lines with arrows represent the electric flux (i.e. the Wilson lines), which 
has to escape the box in a symmetric way through the two opposite planes because of the 
boundary conditions. 


symmetric with respect to 'l'c(ic) but, as discussed in section 6, it might be more practical 
to use in numerical simulations, especially in the framework of compact QEDq. 

Another choice that might look more convenient because of its explicit 0(4) covariance is 
given by 

J^{x) = d^^{x) , df,d^^{x) = 5'^{x) , (3.10) 

where 4>(x) is anti-periodic in space and has appropriate boundary conditions in time. With 
this choice the operator Tj(x) can be written as 

$^(x) = e-nSd^ydpA,{y)Hy-^) . (3.11) 

In Landau gauge we get (x) = V'(x), and the operator ^i{x) is the unique gauge-invariant 
extension of the operator ■!/)(x) defined in Landau gauge. Even though the Landau and other 
covariant gauges are often used in perturbative calculations, notice that the operator T£(x) 
is non-local in time and interferes with the dynamics by effectively generating a time- 
dependent contribution to the Hamiltonian. One can show that this contribution vanishes 
at large time separations, and therefore the same masses will be obtained, but in practical 
situations the asymptotic behavior could be reached very slowly. These complications can 
be avoided in the first place by sticking to a gauge-invariant formalism with the local-in-time 
interpolating operators introduced before. 
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4 Flavour symmetry in QCD+QEDc 


QCD is coupled to QED in the standard way 

{ Nf 'I 

+ | , (4-1) 

where the chromo-magnetic held strength and the covariant derivative are 
Gfj,u{x) = dfj,B^{x) - dyB^{x) - i[B^{x), By{x)] , 

Dl=*d^- iqfA^ - iB^ , (4.2) 

and Bn(x) denotes the colour gauge held. B^{x) is dehned to be a traceless hermitian 
3x3 matrix. Up-type and down-type quarks have electric charge qj = 2/3 and qf = —1/3 
respectively. Since quark helds obey C* boundary conditions, the colour gauge held must 
obey C* boundary conditions as well in order to ensure periodicity of the action density, 

Bp{x + Li) = -Bpix)* . (4.3) 

Let us now focus on the violation of havour and electric-charge conservation, since they are 
substantially different from the case of QEDc alone. 

Since the elementary charge is 1/3, from the discussion in section 2 it might seem that 
processes with a AQ = LLlj?) violation are allowed by the boundary conditions. However 
fractional charges are conhned in hadrons which have integer electric charge. If the box size 
is large enough only colourless particles can travel around the torus, implying that charge 
violation can be produced only in multiples of AQ = ±2. Consequently a proton state can 
mix with an antiproton state, or with a pTr+TT"*" state. 

One might wonder whether C* boundary conditions can induce a spurious mixing of the 
proton with some lighter state. This issue is surely relevant if one wants to extract the 
proton properties from the long-distance behaviour of two-point functions from lattice sim¬ 
ulations. It is also not entirely trivial, considering that C* boundary conditions produce a 
violation of the baryon-number conservation. When a hadron travels around the torus its 
baryon number changes sign, which in turn implies that baryon-number violation can be 
produced only in multiples of 2. A proton state cannot mix with states with zero baryon 
number, i.e. with lighter states. 


Both charge and baryon number are linear combinations of the individual species numbers, 
which we refer to as flavour numbers. 




(4.4) 
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Figure 3: Schematic representation of a possible process responsible for the mixing. 
The process goes throngh a uu pair creation. The colonrless K~ = su travels aronnd the 
torns and tnrns into a K~^ = su. Finally an ss pair annihilates. 


Since each flavonr-nnmber conservation law is violated by the C* bonndary conditions, one 
might wonder for instance whether a pion state can mix with a kaon state. This is not the 
case, as individnal flavonr conservation can be violated again only in mnltiples of two, 

AFf = 0 mod 2 . (4.5) 

Also notice that AB being a mnltiple of 2 implies that total-flavonr F = '^j-Ff violation 
is prodnced only in mnltiples of six, 

AF = 0 mod 6 . (4.6) 

For instance, if only strangeness conservation is violated in a given process, this violation 
mnst be prodnced in mnltiples of 6. If strangeness violation amonnts to a mnltiple of 2 
which is not a mnltiple of 6, then it mnst be accompanied by violation in the conservation 
of some other flavonr. For example the 14“ = sss will mix, via a K~ = su traveling aronnd 
the torns, with + 27 where S'*' = suu and with other two particle states like This 

process has AFg = —2 and AF^ = +2. In particnlar this implies that the 14“ mass cannot 
be extracted from the long-distance behavionr of a two-point fnnction at hnite volnme. In 
order to extract the 14“ mass one has to take the inhnite-volnme limit of the two-point 
fnnction (or effective mass) hrst, and then extract the long-distance behavionr. Similarly 
the H“ = ssd mixes with the p = uud via a K~ = su traveling aronnd the torns (see 
hgnre 3). This process has again AFg = —2 and AF^ = -\-2. 

In QCDc alone, flavonr violation is an exponentially-snppressed effect in the size of the 
box, like any other hnite volnme correction. Adding electromagnetic interactions make 
hnite volnme corrections generically inverse powers of L, dne to the massless photon. The 
detailed analysis of havonr violating process in QCD-|-QEDc reqnires to keep track of the 
havonr nnmbers in the process. This analysis, in the framework of an effective held theory 
of hadrons, is carried ont in detail in the appendix A, bnt the main resnlts that we prove in 
this appendix can be easily explained. Flavonr violating process in QCD-|-QEDc cannot be 
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mediated by the photon. A particle with the same flavour numbers that are violated must 
travel around the torus, and since only massive particles carry flavour in QCD+QEDc, 
these effects are exponentially suppressed. 


For example, in the case of the already-mentioned mixing between the and the proton, 
the one-loop diagram of figure 3 is of order exp(—But the general case is much 
more complicated, since the can also mix with the proton and an arbitrary number 
of photons, or with a neutron-vr"’' state. As it is proved in appendix A, flavour violating 
process in this case are suppressed by a factor exp(—/rL) with 


= 


M; 




/M|_ - Ml, + 

V 2Ms- 


11/2 


(4.7) 


Note that this effects are generically very suppressed, since the corresponding coefficient in 
the two-point function is proportional to the square of the transition amplitude, i.e. to 
exp(—2/iL) ~ 0(10“^*^). A similar analysis for the case of the mixing of the 0“ results in 
an amplitude suppressed by a factor 0(10“®). 

We close this section by remarking that the renormalization of QCD-|-QEDc is not affected 
by electric charge and flavour breaking effects discussed in this section. Indeed these are 
induced by the boundary conditions and locality guarantees that the ultraviolet structure 
of the theory is independent of them. This applies both to the couplings of the Lagrangian 
and to the renormalization constants and mixing coefficients of any composite operator. 


5 Finite-volume effects on the masses of charged hadrons 


The finite-volume corrections to the mass of a stable hadron of non-vanishing charge g, 
which is valid only at first order in and up to corrections in the size of the box that fall 
off faster than any power, can be written as 

Am(L) f gHW , _^ W tTlM t m 4-Mil j- 

m dvr I 2mL 7r(mL)2 AttuiL^ ^ J ’ 

(5.1) 

where m is the particle mass in infinite volume, and m{L) = m + Am(L) is the particle 
mass in finite volume. Typical examples of stable hadrons to which this formula applies are 
the proton, the neutron, the charged pions, the charged kaons, D and B mesons. 

The derivation of eq. (5.1) is given in appendix B. Here we discuss the structure of eq. (5.1) 
that is in fact very simple. 7^ is the ^-th derivative with respect to of the (infinite- 
volume) forward Compton amplitude for the scattering of a photon with energy |k| on the 
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1C* 

2 C* 

3C* 

e(l) 

-0.77438614142 

-1.4803898065 

-1.7475645946 

^(2) 

-0.30138022444 

-1.8300453641 

-2.5193561521 

^(4) 

0.68922257439 

-2.1568872986 

-3.8631638072 


Table 1: Values of the first three coefficients ^(s) in the case of C* boundary conditions 
in 1, 2 or 3 spatial directions and periodic boundary conditions in the others (columns 2,3 
and 4 respectively). 


charged hadron at rest, in the limit k —)• 0. The boundary conditions enter only in the 
dehnition of the generalised zeta function 


f(o) = E 

n^O 


Ini* 


(5.2) 


This formula is valid for real s > 3, while the values s = 1 and 2 are obtained by analytic 
continuation. An explicit representation of the ^(s) coefficients, which is valid for all values 
of s we are interested in, is given in eq. (B.34). The values of the first three coefficients ^(s) 
are given in table 1 in the case of C* boundary conditions in 1, 2 or 3 spatial directions. 


The 1/L and 1/L^ terms are universal, i.e. they depend only on the mass and charge of the 
hadron, and not on its spin and internal structure. The dependence upon spin and internal 
structure is encoded in the coefficients 7} and is suppressed with respect to the universal 
part, as it contributes at 0(1/L^). No inverse odd power of L appears in the expansion, 
other than the leading 1/L point-like contribution. 


A formula very similar to eq. (5.1) has been derived in refs. [2, 15, 16] in the case of 
QEDl, i.e. the theory with the quenched spatial zero-modes of the electromagnetic held. 
According to refs. [2, 15, 16] the 1/L and 1/L^ terms are universal also in QEDl, while 
spin and structure-dependent terms contribute at 0(1/L^). 


It is remarkable that, because of the locality of QEDc, spin and structure-dependent con¬ 
tributions are much more suppressed with respect to QEDl. Moreover, also the universal 
1 /L and 1/L^ contributions are considerably smaller in QEDc with respect to QEDl. This 
is shown in hgure 4 where the results of refs. [2, 15] are compared with the 1/L and 1/L^ 
terms of eq. (5.1). 
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Figure 4: Leading finite-volume corrections to the mass of a stable particle of charge qe 
in QEDc and QEDl- The plot shows the universal 1/L and 1/L^ contributions. Spin and 
structure-dependent contributions are O (l/L^) in QEDl, and O (l/L^) in QEDc- 


6 Lattice formulation 


In the context of lattice non-compact QED, the implementation of C* boundary conditions 
and of the proposed interpolating operators is straightforward. One can extract the leading 
order O(e^) electromagnetic contributions either with the techniques described in [3], or 
by doing a QED dynamical simulation as suggested in refs. [2-7]. In the non-compact 
formulation a way to damp the longitudinal modes of the gauge field is needed, gauge¬ 
fixing being the most common choice. Here we will focus on the compact formulation of 
the theory in a manifest gauge-invariant way. 

In the compact formulation on the lattice, the gauge field is replaced by the link variable 
U{x,ij) which lives in the gauge group U(l) and satisfies the boundary conditions 

U{x + Lk,p) = U{x,p)* , (6.1) 

along the spatial directions, and generic boundary conditions (i.e. periodic, SF, open, open- 
SF) along the temporal direction. For the moment we focus on the simpler case of QEDc 
coupled to a single fermion field with unitary charge. The generalization to QCD-|-QEDc 
will be discussed at the end of this section. 

We want to argue now that, in order to be able to discretize the interpolating operators 
proposed in section 3 in a completely gauge-invariant fashion, we need a rather unconven¬ 
tional action for compact QEDc. Notice that in the standard formulation of compact QED, 
the perturbative series is generated by identifying U{x,p) = and by expanding in 

powers of the gauge field. The discretization of the interpolating operator (3.9) would need 
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to take the square root of the link variable. This operation is not gauge covariant, and 
should be avoided while aiming at a completely gauge-invariant formulation. The root of 
this complication lies in the fact that, because of the boundary conditions, the electric flux 
generated by a single charge must escape the box in a symmetric way through the = 0 
and Xk = L planes. The dynamical unit charge generated by the interpolating operator of 
eq. (3.9) is located in x and sees effectively two image half charges located in x + Lk and 

X - Lk- 


As it will be clear by the end of this section, this issue is completely removed by choosing 
the following action for compact QEDc with a single matter field. 


X flU 

Sm = '^'ip{x)D[U‘^]'ij;{x) . ( 6 . 2 ) 

X 

The plaquette is defined as usual, 

P(x, = U{x, fj,)U{x + fi, u)U{x + V, fi)~^U{x, v)~^ , (6.3) 

while the Wilson-Dirac operator has an unconventional coupling to the gauge field, 

D[u^] = "^ + J E - v;;[c/^]v^[c/2]} , 

^ ti=0 

V^[i7^]V’(x) = U (x, + fi) - ip{x) , 

V))[C/^]V'(a;) = V’(a^) -U{x-fi, - fi) . (6.4) 

Any other discretization of the Dirac operator (preserving charge conjugation) can be em¬ 
ployed as well. 


The proposed action is invariant under local gauge transformations of the form 

U (x, y —)■ A{x)U (x, /r)A(x + fL)~^ , 

V’(x) —)• A{x)‘^'tp{x) , 

Vj(x) —)• '0(x)A(x)“^ , (6.5) 

where A(x) G U(l) satisfies boundary conditions 

A(x -I- Lk) = A(x)* . (6.6) 

The action possesses also a Z| center symmetry. For each direction fi one can flip the sign 
of all link variables in the direction /r on the three-dimensional slice defined by x^ = 0 
without changing the value of the action. Before discussing the interpolating operators, we 
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want to show that the action (6.2) is perturbatively equivalent to the usual QED action in 
the continuum limit. 


In order to set up a perturbative expansion, we need to identify the minima of the action at 
O(e^). These are given by all configurations with P{x,^,v) = 1. In appendix C we show 
that there is a discrete set of gauge-inequivalent minima labeled by the elements of the set 

n = {(zo, Zi, Z 2 , 1) \ Zq = zf = Z 2 = 1} ■ (6.7) 

Given a minimum of the action at O(e^), it is always possible to find a vector z £ Q such 
that the chosen minimum is gauge-equivalent to the following gauge field 


Uzix,fJ,) 



a Xf, = - 1 , 

otherwise . 


( 6 . 8 ) 


Because of center symmetry one might expect also minima with z^ = —1. However each 
minimum with Z 3 = —1 is gauge-equivalent to some minimum with Z 3 = 1 (this is a 
byproduct of the construction given in appendix C). 


The perturbative expansion around the minimum Uz{x,^) is set up by defining 

U{x,n) = C72 (x,, (6.9) 

and by adding a gauge-fixing term Sgf to the action, which we will not do explicitly. We only 
observe that Sgf is a function of the fluctuation A^{x) only, and not of the classical vacuum 
Uz{x,fj,). Given a generic functional of the fields, the perturbative expansion to 

some order 0 {e"') of its expectation value is given by 

= e-5[e2^.b.b]-5gf[A] ^ ^ 

( 6 . 10 ) 

where we have used that the action is center-invariant and therefore it does not depend on 
Uz{x,fj,) once the substitution (6.9) is used. The normalization Z is given by 

Z = 8 J VAVtpV^p {1 + C)(e’")} . (6.11) 

If the observable J- is charged under center symmetry, its expectation value vanishes. On 
the other hand, center-invariant observables get mapped naturally into corresponding ob¬ 
servables in the non-compact setup. In fact, if T is invariant under center symmetry, then 
the 2 ; dependence drops out of the path integral and the standard perturbation expansion 
about ?7 = 1 is recovered, 

= ^ J VAVi^V^ + ©(e*") , (6.12) 
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provided that S[e ^^is the standard QED action up to irrelevant operators. 


This can be verified by replacing the definition (6.9) into the action and by expanding in 
powers of the fields. The 1/2 factor in the exponent of (6.9) combines with the uncon¬ 
ventional normalization of the gauge action in eqs. (6.2) in such a way that the canonical 
normalization of the gauge field is restored, 

P(x, ^l,v) = l + ^F^u{x) - + ... , 

<57 = ^ ^ ^[1 - Fix, ^ X] X] F^u{x) + irrelevant operators . (6.13) 

X fiv X flU 

Also the same 1/2 factor in the exponent of (6.9) combines with the second power of the 
link variable in the Dirac operator (6.4) in such a way that the correct coupling of the 
electron to the gauge field is restored, 


U (x, = 1 + iAfj,{x) + 


Sm = Yl '^(®) I > 


9u + d* 


+ iA^{x) 


+ m> 'ip{x) + irrelevant operators . (6.14) 


The elementary charge (charge is quantized in compact QED) interacts with the gauge 
field with strength 1/2. However the dynamical fermion has an electric charge that is twice 
the elementary charge, which generates a coupling of strength 1 to the gauge field. This 
structure is also reflected by the gauge transformations (6.5). 


In the proposed setup, thanks to the identification (6.9), the interpolating operator (3.9) 
can be discretized in a straightforward fashion, 

— 1 L—x/c —1 

'^' 3 ( 3 :) = U{x + sk,k)~^'ijj{x) U{x + sk,k). (6.15) 

s=—Xk s=0 

Notice that the the above operator is charged under center symmetry. However in practice 
only the product 4's(x)Ts(y) is relevant, which is center invariant. 


For completeness we present also a possible discretization of the operator (3.7). We intro¬ 
duce the field 

Alix) = A-^dlFf,^ix) , (6.16) 

where A = 9^5^ is the three-dimensional discrete Laplace operator defined with anti- 
periodic boundary conditions, and Fpa{x) is some discretization of the field tensor (e.g. the 
clover plaquette). It is straightforward to verify that A))(x) satisfies the discrete Coulomb 
constraint d^A^{x) = 0. In the continuum limit A'^(x) is nothing but the gauge field in 
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(6.17) 


Coulomb gauge. ^ The operator (3.7) can be discretized by using the relation 

Tc(a:) = ^'5(2:) , 

which is exact in the continuum limit, and easily verified in Coulomb gauge. In fact in 
Coulomb gauge and in the continuum eq. (6.17) is completely equivalent to eq. (3.9), given 
the relations A^{x) = A'^{x) and Tc(x) = 'ijj{x). 


The generalization of the proposed strategy to the case of compact QCD+QEDc is straight¬ 
forward. We need to introduce the link variables V{x,^) G SU(3) for the colour field with 
the boundary conditions 

V{x + Lk,p) = V{x,py , (6.18) 

and the corresponding plaquette: 

Q{x, p,,v) = V{x, p)V{x + (1, i')V{x + u, p)~^V{x, iy)~^ . (6.19) 

For sake of simplicity we choose the standard Wilson action for the colour field. The photon 
action requires a further rescaling, since quarks have fractional charge. 


S — Sq Sm ; 


Sg = ^ 

g 
18 

e 


sEE tr [1 - Q{x,p,u)] , 


X flV 

X fll/ 

/ ^ 


( 6 . 20 ) 


Moreover the Dirac operator has to implement the correct coupling of up-type {qf = 2/3) 
and down-type {qf = —1/3) quarks to the electromagnetic field, 

1 ^ 

Df[U, V]=mf + -Y, {7m(V* - V*yU^'ifV]V g[U^'ifV]} . 

^ ti =0 

( 6 . 21 ) 


We remind that QCD with non-degenerate Wilson-Dirac quarks (with or without QED) 
has a mild sign problem, i.e. the fermionic determinant is positive in the continuum limit 

® In the continuum limit: 

dkAl{x) = A~^dkdjFkj{x) = 0 , 

A1{x) = A~^dk{dkAf,{x) - d^Ak{x)} = A^{x) - {A~^dkAk{x)] , 
i.e. T“(a;) is gauge-equivalent to T^(a:) and satisfies the Coulomb-gauge contraint. 
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but can be negative because of lattice artefacts. In appendix D we show that C* boundary 
conditions do not make this sign problem worse. 


7 Conclusions 


A local solution to the problem of electrically charged particles in a finite volume was 
proposed in [18-20, 26], and it is based on C* boundary conditions for all fields along one 
or more spatial directions. Because of the boundary conditions Gauss’s law does not prevent 
the propagation of charged particles on a finite volume (as opposed to the case of periodic 
boundary conditions). We have analyzed in detail the properties of QED in isolation and of 
QED coupled to QCD with C* boundary conditions (QEDc and QCD+QEDc respectively), 
and we have discussed how this setup can be used in spectroscopy calculations. 

We have devoted part of the paper to construct interpolating operators that have the 
quantum numbers of charged particles and that are also invariant under local gauge trans¬ 
formations. These can be used to probe the physical sector of the Hilbert space of the 
theory with non-perturbative accuracy without having to rely on gauge-fixing at interme¬ 
diate stages of calculation. To this end we have discussed the details of the implementation 
of the proposed interpolating operators in the compact lattice formulation of the theory. 

We have discussed the symmetries of QEDc and QCD-|-QEDc in depth. In particular we 
signal that C* boundary conditions violate flavour and electric-charge conservation partially, 
in such a way that this does not represent a limitation to the use of C* boundary conditions 
in most of the relevant applications. Even though finite-volume effects vanish generally like 
some inverse power of the box size because of the photon, we have shown that flavour and 
electric-charge violations are exponentially suppressed in the box size. 

We have calculated the finite-volume corrections to the masses of charged particles with 
C* boundary conditions at O(aem)- We have shown that the leading 1/L and 1/L^ finite- 
volume corrections are universal, i.e. they depend on neither spin nor internal structure. 
Similar results have been previously obtained in the non-local formulation QEDl. When 
compared with these previous result, the finite-volume corrections of QEDc are found 
to be significantly smaller. In particular, the non-universal spin and structure-dependent 
corrections are 0(1/L^) in QEDl and 0(1/L^) in QEDc. We have also shown that these 
non-universal terms are related with physical quantities, namely the derivatives of the 
forward Compton scattering amplitudes. 
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A Exponential suppression of flavour mixing 


Let H(x) be some interpolating operator for some fixed spin-component of the negatively 
charged = ssd and H_|_(x) its C-even component. We consider the finite-volume 
Minkowskian retarded two-point function at zero momentum, its spectral decomposition 
and the dispersion relation: 

C{E-,L) =t [ d\ 0(xo)e*®"°(H+(x)tH+(O)) = H df, , (A.l) 

J]Rxi,3 Jo fi-E-ie 

p{E-,L) =-lmC{E;L) . (A.2) 

TT 

in infinite volume the spectral density vanishes for E < M^-. In a finite box, because 
of C* boundary conditions, the lowest state contributing to the two-point function is a 
proton state (via a strangeness-violating process) and the spectral density vanishes only for 
E < Mp{L). We want to show that the spectral density vanishes exponentially with the 
volume for energies lower than M=-. 


More precisely we choose some energy E < M'^- and a smooth test function (pEip) which 
vanishes for p, > E. Then we want to show that, in the L —)• oo limit 


In 



dp p{p; L)(I)e{p) < —2LM{E) + C>(lnL) , 


(A.3) 


where the mass that controls the exponential decay is a decreasing function of E < M=- 
and 


M{E) > 




M| - Ml, 


+ M‘^ 




2M=- 


11/2 


(A.4) 


Notice that the spectral density at finite volume is a sum of delta functions localised on 
the eigenvalues of the Hamiltonian in the given channel,® which is the reason why we need 
to consider a test function in order to write a precise statement. Before proceeding we 

®At fixed order in perturbation theory, the spectral density is a sum of delta functions and their deriva¬ 
tives, localized on the eigenvalues of the free Hamiltonian 
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comment on the fact that the analysis presented in this appendix can be easily extended 
to other channels (e.g. to the mixing of the with lightest states). 

We assume that the leading finite-volume corrections in the two-point function C{E\L) are 
described by some arbitrarily-complicated Lagrangian field theory with small couplings, 
which effectively describe the dynamics of hadrons and photons at large distance in the 
framework of a perturbative expansion (after gauge-fixing for the photon). In order to 
avoid IR divergences at any stage of our calculation, we assume that the photon is massive. 
We will find that the mass A4(E) does not depend on the mass of the photon. Each 
infinite-volume stable particle is described by an elementary field in the effective Lagrangian. 
Because of locality the finite-volume theory is described by the infinite-volume Lagrangian 
density. In particular vertices conserve flavour. Fields are assumed to have definite flavour 
numbers, and fields with all flavour numbers equal to zero are assumed to have definite 
C-parity. 


We think of the two-point function C{E]L) order by order in perturbation theory as a 
sum of Feynman diagrams. According to Cutkosky’s rules, a Feynman diagram contributes 
to the spectral density p{E]L) via the dispersion relation (A.2) only if a cut between its 
external vertices exists such that the sum of the masses of the cut propagators is smaller 
than E. Notice that all states propagating between the two external vertices must have 
(—1)^ = —1 and ( — 1)'^ = 1 where B is the baryon number and S is the strangeness number. 
It is easy to check that (at the physical masses) states with such quantum numbers and 
with energy lower that M=- can contain no strange particle. Therefore a Feynman diagram 
for the two-point function contributes to the spectral density for E < M=- only if a cut 
exists between the external vertices such that no strange particle propagates through the 
cut. The set of these diagrams is denoted by Eg (see fig. 5 for an example). The spectral 
density for E < M=- can be represented as 


p{E;L) 


— im 



0(xo(ei))e*''"°«^)Ff (x)U(^„)=o , 

(A.5) 


where and ^-re the two external vertices of the diagram and V{Q) is the set of all 
vertices. E^(x) is a function of the coordinates of all vertices of the graph, and it is given 
by a product of propagators in coordinate space, various derivatives of propagators and 
numerical coefficients. We will refer to the function F^{x) as Feynman integrand. 


Structure of the Feynman integrand. We denote by L{Q) the set of lines of the 
diagram Q. Each line i originates from the vertex i{i) and terminates in the vertex f{i) 
(line orientation is chosen arbitrarily). Each line i corresponds to the Wick contraction 
of two fields located in i{i) and f{i). Let E{£,i{i)) and F(i,f(i)) be the vectors that 
contain all flavour numbers of these two fields, F = {U, D, S,...). It is convenient to define 
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Figure 5: Example of a Feynman diagram in Vg. This is a contribution to the 
(H+(x)^H+(0)) two-point function. Propagators with a cross are of the flavour-violating 
type, and they exist only because of the C*-boundary conditions. Either the proton or the 
pair n -|- tt'*' (in both cases with an arbitrary number of photons) can go on-shell with an 
energy lower than M^-. 


F{i, a) = 0 if the line i is not attached to the vertex a. Since a field can be Wick-contracted 
either with itself or with its charge conjugate, each line is uniquely associated to a mass 
m£. The line I corresponds to a propagator 

A^^{5ex-,me,L) = E Cf (n^) AuiSexo, 5ix + Lrit, me) , (A.6) 

where AM{x,m) is the infinite-volume Minkowkian propagator, 5ex = x{f{i)) — x{i{i)), 
ne = (O,!!^), and C'^{'n.e) is a function defined by one of the following possibilities: 

Type 1. Line (. arises from the Wick contraction of two flavourless C-even fields, 

= F{i,f{i)) = 0 , Cfine) = 1 . (A.7) 

Type 2. Line i arises from the Wick contraction of two flavourless C-odd fields, 

F{e,i{i)) = F{iJ{i)) = 0 , cfine) = . (A.8) 

Type 3. Line i arises from the Wick contraction of two flavourful fields with opposite 
flavour, 

Fii,iii)) = -FiiJii)) / 0 , Cfine) = . (A.9) 

Type 4. Line i arises from the Wick contraction of two flavourful fields with same flavour, 
Fi£,iie)) = FiiJii)) / 0 , Cfine) = . (A.IO) 

Notice that conservation of flavour at the internal vertices is expressed by the equation 

^F(£,a) = 0 , for aGF(g)-{eo,6} • (A.ll) 

e 
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(A.12) 


The function (x) has the following structure 

n 

F^{x,n)=Y^ Cf{n() AM{Sex + Lni;me) , (A.13) 

eeL(g) 

where is a differential operator acting an all the coordinates x{a) (in fact derivatives 
might be inserted in between propagators), and includes also couplings and combinatorial 
factors. The detailed structure of is of no interest for the current discussion. 


Boundary conditions for the Feynman integrand. C*-boundary conditions for the 
helds imply some peculiar boundary conditions for the Feynman integrand. Shifting the 
coordinate of a single vertex by 

x{a) —)• x{a) + Li (A.14) 

is equivalent to replacing the operator inserted at the vertex a with its charge-conjugate 
(the propagators have to be modihed accordingly). The diagram obtained by this procedure 
is denoted by CaiG)- Since the interaction Lagrangian and the considered interpolating 
operators are invariant under charge-conjugation, ^ is a diagram contributing to the two- 
point function if and only if CaiG) is a diagram contributing to the two-point function. It 
is also easy to check that the class T>s is closed under the action of Ca- By iterating the 
action of Ca, and by noticing that is the identity we get 

Ff(a:o,x + LA) = Fj;,j;,(a:o,x) , A(o) e Z“ , (A.15) 

c^(e)= n Y"He). (a.i6) 

a&V{g) 

Under the action of Ca, the vertex operator does not change. Propagators of type 1 attached 
to the vertex a (at only one of the endpoints) are invariant, propagators of type 2 flip sign, 
propagators of type 3 are replaced with propagators of type 4 and vice versa. In formulae 
this is equivalent to: 


< 

S 

II 

< 

(A.17) 

= Cfiw - SeX) . 

(A.18) 

The boundary conditions for the function F^{x,n) follow 


F^(xo,x-hLA,n) = F^^g^{xo,:x.,n + 5X) . 

(A.19) 


We can use the above boundary conditions to restrict the sum over the pairs (^,n) and 
simultaneously extend the coordinate integration range to the whole M^, by means of a 
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construction that has been discussed already in [22], We will say that the two pairs {Q, n) 
and {Q', n') are gauge-equivalent if and only if A(a) G with a G V{G) exists such that: 


G' = c^iG) , (A.20) 

= n£ + 6i\ = n£ + A(/(£)) - A(i(^)) . (A.21) 


The field n defined on lines is referred to as gauge field, and the field A defined on vertices is 
referred to as gauge transformation. The set of all possible pairs {G, n) splits into equivalence 
classes denoted by [(^,n)]. As shown in [22], given two equivalent gauge fields the gauge 
transformation that relates the two of them is unique up to a global gauge transformation. 
The sum over the pairs (^,n) can be written as a sum over the equivalence classes and a 
sum over the gauge transformations with A(^o) = 0. The spectral density for E < M=- 
becomes 


p{E;L) 


— im 
vr 

[i5,n)] 

s.t. 



(A.22) 


By using the boundary conditions (A. 19), one can use the sum over the gauge transforma¬ 
tions to reconstruct the integrals over 


p{E-,L) = -lm 

TT 


[(^,n)] 
s.t. Q^i^s 



d^x{a) > X 




(A.23) 


Strangeness flow in Feynman diagrams. We introduce some general definitions. 

Paths. A path P connecting the vertices a b is a set of lines, with the property that a 
sequence a = vi,V 2 , ■. ■ ,vm = ^ of pairwise different vertices exists, together with a labelling 
of all lines in P, such that Vk, r’fe+i are the endpoints of ik- 

Loops and Wilson loops. A loop C is a set of lines, with the property that a sequence 
vi,V 2 , ■ ■ ■, V]\f of pairwise different vertices exists, together with a labelling ii,i 2 , ■ ■. ,iN of 
all lines in C, such that v^, are the endpoints of for k = 1, • • • , A^ — 1, and vjsf, vi 
are endpoints of ii\f. The Wilson loop associated to C is defined by 

W{C,n) = (A.24) 

l£C 

Notice that a Wilson loop is gauge invariant. 
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Trees and axial gauge. A tree T in a connected graph is a maximal set of lines that contains 
no loops. If 7~ is a tree of G, for any pair of vertices a ^ b of G, there is a unique path 
P C 7~ connecting a and b. Any connected graph with more than one vertex has at least 
one non-empty tree. Given a tree T, a gauge field n is said to be in axial gauge with respect 
to T if = 0 for any i € P. It is easy to show that any gauge field is gauge-equivalent to 
some gauge field in axial gauge with respect to T. 

We give a closer look at the flavour structure of diagrams in focusing in particular on 
strangeness. A line i is said to be strange if S{i,i{i)) ^ 0, and strangeless otherwise. For 
any diagram G & Ps, define the subdiagram Gs by taking only the strange lines in G 
and vertices that are attached to these lines. Clearly the vertices and ■?! corresponding 
to the interpolating operators belong to Gs- Because of the defining property of 2?^, there 
is no path in Gs connecting and ^i. We define as the connected component of Gs 
containing ^o- By specializing eq. (A.11) to strangeness, and by observing that strangeness 
cannot flow outside of Gs,o, one obtains the following: 

Proposition A.l. Strangeness is conserved within Gs,o o-t M vertices ofGs,o except i.e. 
^ S'(2,a) = 0, for a eV{Gs,o) - ^0 ■ (A.25) 

ieL{gs,o) 


Lemma A.2. For any pair (^,n) such that G ^ Ps F^{x,n) is not identically zero, 
then a loop C exists in Gs,o such that |VF(C', n)| > 1. 


Proof. We consider a tree 7~ in Gs,o, and we assume without loss of generality that (^,n) 
is in axial gauge with respect to the tree P. Notice that all lines in ^^,0 are flavourful, 
therefore they can be either of type 3 (flavour-preserving) or of type 4 (flavour-violating). 
If 2 G T then = 0 by definition of axial gauge. In this case i must be flavour-preserving, 
otherwise the propagator would contribute with a (5(n^)q = 0 factor to F^(x,'n). 


If all lines in Gs,o were flavour-preserving, by conservation of strangeness at the internal 
vertices, eq. (A.25), we would have 


0 = 


E 


E 


a6V'(5s,o)—Co t£L(gsft} 


^ ^ S(£,a)- 

<isy(Ss,o) ^SL(Gs,o) i&L(Qsfi) 


= Y. m,m+s(flf(i))]- Y s(iPo) = - Y *^(^’^ 0 ), 

t&pGsfi) eeL{g,,Q) teL(g.,o) 


(A.26) 


which is in contradiction with the fact that the interpolating operator has strangeness equal 
to ±2. Therefore at least a flavour-violating line i G L(Gs,o)—P exists. Since the propagator 
of a flavour-violating line comes with a factor then necessarily 0. 
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Consider the path P in T that connects i{tj to f{i), then it is straightforward to prove 
that C = PL) {£} is a loop in Qs,o 

|VF(C,n)| = |n^-| /O . (A.27) 

The thesis follows from the fact that the components of are integers. 


□ 


Time-ordering and Euclidean kernel. It is convenient to separate different time- 
orderings in integral (A.23). Formally the set T{Q) of all time-orderings of G is defined 
as the set of permutations of V{G), i.e. the set of all bijective functions between V{G) 
and {1,2,...,|E(^)|} where |E(^)| is the number of vertices of Q. Given a time-ordering 
r G T{Q), and two vertices a,b € V{G) such that r(o) < r(6), we will say that a is before b 
and b is after a (with respect to r). We define the time-ordering function associated to r 
as 


Ot{x) = 6{xo{b) - xo{a)) , (A.28) 

a,bevig) 

s.t. T(a)<r(fo) 

and we insert the identity in the integral (A.23) in the form of 

1 = Orix) . (A.29) 

TeT(e) 

Not all time-orderings contribute to the spectral density for E < M=-. States with no 
strange particles need to be able to propagate at some intermediate time between the two 
external vertices. Therefore only time-orderings such that is after any vertex in Gs,o 
contribute. We define Ts{Q) the set of such time-orderings. Eq. (A.23) becomes 

/9(E;L) = ilm Y M 11 ld^x{a)\x 

^ [(g,n)] T&Ts(g) J 

S.t. Q^i^s 

X . (A.30) 


We choose a diagram and a time-ordering contributing to the previous formula. Without 
loss of generality we can assume that each line is oriented in such a way that its final point 
is not before its initial point. Let vf be the latest vertex in Qsfl, i.e. the vertex in Qsfl with 
the property that any other vertex in is before it. Let v'!^ be the vertex right after vf, 
i.e. the vertex defined by = T{vf) + 1. 

We construct the subdiagram Qf (resp. QZ.) by taking all the vertices in Q not before v'f, 
(resp. not after vZf) and all lines in Q connecting any pair of these vertices. Let Lg be 
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the set of lines with one endpoint in Q'[_ and one in Q'!^. This decomposition induces a 
factorization of the integrand in eq. (A.30). The time-ordering function factorizes as 


Orix) = 0{xo{v\) - Xo{v'[_))er+{x)Or_{x) , 
Ot±{x)= e{xo{h) - xo{a)) , 


(A.31) 

(A.32) 


a,b&V{gi) 

s.t. r(a)<r(b) 


where r+ (resp. r_) is the time-ordering restricted to Q'!^ (resp. QZ.). The Feynman 
integrand factorizes as 

Fg{x,^) = (x,n) C'f(n£) A+-h Ln^; m^) , (A.33) 

where /x is some collective Lorentz index (in the following we will omit the sum over ^). No¬ 
tice that Flfgr {x, n) (resp. FjfgT {x, n)) contain differential operators acting on coordinates 
of the final (resp. initial) points of the lines in Lg. A^{x,m) is the retarded Minkowskian 
propagator which we can conveniently write in time-momentum representation 


A+(x;m) = e{xo) j , with E = + . 


(A.34) 


We plug the factorizations (A.32) and (A.33) and the explicit representation of the retarded 
propagator into eq. (A.30), we substitute xo{a) —)• xo(a) + xo(x-) for any a G V{Q+), and 
we use invariance under translations of F)(^^(x,n). The integrals over the coordinates 
factorizes over the two subgraphs: 


= E i 11 


d?p(. 


[(e?,n)] t^Ts{G) [i&L- ^ ^ \ 

s.t. 

X i‘^T^f^^iE£eL^Pi)K~g^^{E,p,n)R^g^^{E,p,ii) . 

(A.35) 

The explicit expression for R^g t{E, p, n) is given for sake of completeness, but we will not 
need it in the current discussion 


Flg,r(F,P,n) = 


OteLj Ff{ni) 
SfeLj Fi E ze J 


TT [ d'^x{a)\ er+{x)x 

ZTZt \ J** 


X F/^e.(x,n)5^(x(z;;))e*®^°(«i) e-*Ez^o(/W)g*Pzx{/(£)) ^ 


iGLZ 


The explicit expression for K^g ^{E,p,n) is given after Wick rotation of the coordinate 
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integrals xo(a) —)• —ixo(a): 





Or. {x)Fl^^gr_ (a;, n) x 


^ gExo(i;I)g-E£6L5 ■E£[xo(i;I)-xo(iW)]^-iE^6L5 P<.[x(i(£))+Ln<.] 


(A.37) 


where = y rnj + p|, Fqt {x, n) is constructed as {x, n) except that the Minkowskian 
propagators and vertices are replaced by the Euclidean ones, and is some constant 

phase factor. The Wick rotation is allowed for Kg ^{E, p, n) because all states propagating 
in between the interpolating operator and the latest vertex wl have energy that is higher 
than E by construction. We will refer to the function Kg ^{E,p,n) as Euclidean kernel 
associated to the diagram G € T>s and the time-ordering r G Ts{Q). 


We highlight the following observation, which follows from the construction of the Euclidean 
kernel. 

Observation A. 3. The latest vertex v'L of Gs,o is olso the latest vertex of Q'L- Since Q'L 
contains all vertices of Q not after v'L, and all lines of Q with both endpoints not after v'L, 
then Qsfi is a subdiagram of QZ. ■ 


Large-volume behaviour of the Euclidean kernel. We will see that the Euclidean 
kernel provides the exponential suppression in the inhnite-volume limit. We will establish 
this result in a few steps. 

Theorem A. 4. For large L, we have 


= -L£g^r{E,p,K} + 0{\aL) , 
where the function Sg^riE, p,n) is given by 
^e,r(£^,P,n) = 


= min l-Exo{vZ) + '^ E(,[xq{vZ) - xo{i{t))]+ ^ m(,\5iX + n(\ 

( t&Li teL(ei) 

F>g,T = {( 2 :(a) )a&v(gL) I x{a) G = 1 , x(^o) = 0} . 


(A.38) 


(A.39) 

(A.40) 


Proof. We use the heat-kernel representation for the Euclidean propagator 
/■oo 

AE{x;m) = ds — —^ . (A.41) 

Jo (4vrs)^ 

When plugging this in eq. (A.37), we need to introduce a Schwinger parameter si for each 
line i G L{QZ). By substituting si —)• siLj^mi and x(a) —)• Lx{a) we get the general form 
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for the Euclidean kernel 


/*00 _ p 

= n / n / d‘^x{a)6{xo{Co))9r-{x)x 

i.&L(gz)'^^ aey(ei)-& 

ja Pfijs, X, n) ^-L\x{s,x,E,p,n)-zJ2eeLj: Pdx(*W)+nd} 

^ ^ r\( \ ^ ’ 

Q\s) 


(A.42) 


where -P(s, x, n) and Q{s) are polynomials that come from the derivatives in the vertex 
operator and explicit powers of s in the heat-kernel representation, and 


A(s,x,E;,p,n) = 

=-Exq{vZ_)E e,[xo{v'[_) - XQ{i{^))]+ ^ ^ IS£ +—\6iX + ni\'^\ . 

i&Li i^L{g^) I J 

(A.43) 

The large-L expansion of the integral in eq. (A.42) is given by a saddle-point approximation, 
i.e. by expanding about the minima of X. This yields the asymptotic behaviour (A.38) 
with 


£g^T(E,p,n') = min min Ai(s, x,-E, p, n) . 

’ x&Dg^T s^g[0,co) 

Eq. (A.39) is obtained by performing the trivial minimization over each si. 


(A.44) 


□ 


The statement of theorem A.4 makes sense because Eg^riE^Pj'n) is strictly positive. We 
will show this fact in two steps, and we will provide also a lower bound for this function. 

Theorem A.5. The function Sg^r{E,p, n) defined in eq. (A.39) satisfies 

£g^r{E, p,n) > 0 , for E < M=- . (A.45) 


Proof. Let N be the number of vertices of the diagram Q'H. We assume that the orientation 
of lines in Q'H is chosen in such a way that 

rif{i)) > . (A.46) 


We construct an auxiliary graph TL by starting from an empty graph, and by adding elements 
to it accordingly to a set of rules. 


1. We grow TL by adding all the vertices of Qf. At this stage: 

V{Ti) = V{gf), L{TL) = 0. 


(A.47) 



We equip the graph % with the time-ordering function r inherited from Q'L- 

2. For each mother line i of QZ. such that — T{i{t)) < 1, we add the daughter line 

i to % with the same incidence relations (notice that the endpoints of i are already 
in %). 

3. For each mother line i of GZ such that t(/(£)) — r(i(£)) > 1, we imagine to cut the line 

i at each intermediate timeslice. More precisely we add the n = T{f{i)) — — 1 

intermediate vertices vi,... ,Vn and the n -|- 1 daughter lines £o) • • • to % with the 
following incidence relations: 


i(4) = i{e) , 

/(4) = ^>1 , 


(A.48) 

i{£k) = Vk , 

/(4) = Vk+i , for A; = 1 ,.. 

. ,n - 1 , 

(A.49) 

i{£n) — Vn , 

/(4) = f{£) . 


(A.50) 


We also time-order the extra vertices and we declare that 

T{vk) = T{i{l)) + k . (A.51) 

We associate a mass and a gauge field to each line of T-L. The masses of the 
daughter lines are all equal to the mass of the mother line, while the gauge field of 
the original line is transferred completely only to the first daughter line, i.e. 


M4 = mu , 
Wfe = W<5fc,o • 


(A.52) 

(A.53) 


All quantum numbers are also naturally transferred from the mother line to the daugh¬ 
ter lines. 

4. For each mother line ^ of Lq such that t(vZ) — T{i{i)) > 0, we add a final point at 
the latest timeslice and we imagine to cut the line i at each intermediate timeslice. 
More precisely we add the n = t{vZ) — r(i(£)) vertices vi,... ,Vn and the n daughter 
lines f'o,... (-n-i to % with the following incidence relations: 


i(4) = i{tj , /(4) = vi , 

i{h) = Vk , /(4) = Vk+i , for A: = 1, ..., n - 1 

Analogously to the previous case, we declare: 

T{vk) = r{i{£)) + k , 

= El , 

Wfe = w4,o , 


(A.54) 

(A.55) 


(A.56) 

(A.57) 

(A.58) 


and we transfer the quantum numbers of the mother line to the daughter lines. 

The above construction is schematically represented in fig. 6. The auxiliary diagram Ti is 
useful to rewrite the function £g^T-{E, p, n) defined in eq. (A.39) in a way that its timeslice 
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Figure 6: Exemplification of the construction given in the proof of theorem A.5. We start 
with a diagram whose vertices have been time-ordered (top diagram). Each vertex defines 
a timeslice. The diagram in the bottom is obtained from the top one by inserting a new 
vertex every time a line intersects a timeslice. Vertices in the same timeslice are forced to 
have the same time coordinate. 


structure be manifest. 

We define Vk{%) for A: = 1,..., V as the set that contains all vertices a with r(a) = k, i.e. 

Vkin) = T-^k) . (A.59) 

The sets Vk{T-i) for k = 1,... ,N constitute a partition of 
N 

v{n) = □ Vk{n). (A.60) 

k=l 

We will refer to 14(?A) as the k-th timeslice of %. We define Lk{'H) for fe = 1,..., V — 1 
as the set that contains all lines (. such that i{i) is in the A:-th timeslice and f(i) is in the 
(A:-|- l)-th timeslice. The only lines of Ti that are left out are the looplines i with i(i) = f{i) 
which we collect in the set LioopiT-L), 

N-l 

L{n) = □ Lk{n) u Lioopin). (A.ei) 

k=l 


We change the way in which we associate coordinates to vertices. We associate the time- 
coordinate 4 to the timeslice 14 (?A), and independent spatial coordinates x(a) to each 
vertex a G V{'H). We force all vertices in the same timeslice to have the same time 
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coordinate. Therefore the four-dimensional coordinate vector is mapped into 

x(a) = (t^(a),x(a)) . (A.62) 

Notice that in each timeslice Vk{Ti) there is only one vertex of the original diagram 
therefore tk does coincide with the temporal coordinate of this vertex. 


By using iteratively the property that the shortest path between two points is the straight 
line, in the form of 

min I (4+2 - tk+iY + |x 2 - xip {tk+i - 4)^ + |xi - xopj = 

= V (4+2 - tkY + |X 2 - xo|2 , for 4+2 > 4+1 > 4 , (A.63) 


and the property that the shortest path between a point and a plane is the the straight line 
that is orthogonal to the plane, in the form of 

min (4+1 - 4)2 |xi - xqP = 4+i - 4 , for 4+i > 4 , (A.64) 

Xl 

and by recalling that = N as v'[_ is the latest vertex in one can easily prove that 


‘?g,r(4,p,n) = 


= mm mm 


-E[tN -U(^o)]+ Mi\x{f{i)) - x{i{e))+ ne\ 

eeL{n) 


(A.65) 


We choose quantities associated to the lines of Ti satisfying the following constraints 

A£ = 0 , for looplines i G Lioop{T-L) , 

A^ = 0 , for £ G Lk{T-L) with k < t{^q) , 

0 < A^ < Mi , for any £ , 

E = y] Ai , for any k > r(^o) • (A.66) 

teLfc(W) 

Notice that the last two constraints are in tension with each other. Quantities A^ satisfying 
such constraints exist for E < M=- thanks to the following observation. 


Observation A.6. By construction of Qf, only states with energy not lower than M=- 
propagate at any time between the vertiees vf, i.e. 

yy Mi > M=- , for any k > r(^o) • (A.67) 

£&Lk{H) 


Besides the above constraints, the quantities A^ are largely arbitrary (and we will use this 
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arbitrariness later on). The contribution in eq. (A.65) can be distributed over 

the lines of % yielding 

= 

= min min ^ {-/\t[xo{f{t))-XQ{A))\+Mi\x{f{tj)-x{i{t))+ni\], 

(A.68) 

We can provide a lower bound for each term of the above sum by using the following 
inequality 

-A|,5t| + My/ {Sty + > |z| Vm2 - A2 , (A.69) 

valid for any 2 and 6t, as long as M > A. Therefore we get an estimate for £g^T-{E, p, n) in 
which the time minimization has disappeared 

^e,r(A^,P,n) > mm ^ |x(/(^)) - x(z(£)) + n^l - A^ , (A.70) 

i&L{n) 

and which shows explicitly that 

£g^T{E, p,n) > 0 , for ill < M^- . (A.71) 

□ 


Theorem A.7. The function £g^T{E,p, n) defined in eq. (A.39) satisfies 
£g^r{E, p, n) > Ai {E) , for E < M=- , 
where the function Ai{E) is given by 


M{E) 


'Mk± 

[a 4± - 
Mj*- 


E-M, 


1/2 


E^-Mlo+MA 

2E 


1 1/2 


if 0 < E < Mp 

m2 m2 

ifMp<E<^Aur^ 


< g < Mh- 


(A.72) 


(A. 73) 


Proof. We use the construction in the proof of theorem A. 5. Since each term in the sum in 
the r.h.s. of eq. (A.70) is positive, a looser lower bound on the function £g^r{E,p,n) can 
be provided by restricting the sum to a subset of L(fH). 

We choose a strange loop C in with |1T(C', n)| > 1. Such a loop exists thanks to 
lemma A.2. Since is a subgraph oiQ'L, C is also a loop in Q'L. We construct the loop C 
of % by replacing each mother line in C with the corresponding daughter lines. Since the 
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gauge field of each mother line is equal to the sum of the gauge fields of the corresponding 
daughter lines, it is clear that 


|iy(C,n)| = |l^(C,n)|>l. 


(A. 74) 


By using this particular loop in eq. (A.70), we get the following bound for p, n): 

£g^r{E,Y>,T^) >nAn^ |x(/(£)) -x(f(^)) + n(\yjMj - Aj > 

£ec 

>M(E, A) miny^ |x(/(f')) — x(f(£)) + n^l > 
i^C 

>M{E, A)\W{C, n)| > MiE, A) , (A.75) 

where we have used the triangular inequality iteratively along the loop, and we have defined 

M(E, A) = min JMj - Aj . (A.76) 

iec ^ 

The A^’s are arbitrary quantities satisfying the constraints in eq. (A.66). We can use this 
arbitrariness in order to optimize the above lower bound. 

We separate two cases: either C is a single loopline or C has more than one line, none of 
which being a loopline. In the first case, A^ = 0 for £ G C and 

MiE, A) = Mi > Mk± , (A.77) 

where the last inequality comes from the fact that the is the lightest strange particle. 


We consider now the second possibility, i.e. C has more than one line none of which being 
a loopline. 


N-l 

c c y Lkin) . (A.78) 

k=l 

We define 


Ck = cnLkin) . 

The minimization in eq. (A.76) can be performed in two steps, according to 
MiE, A) = min min \ M? — A? . 

^ s.t. i^Ck ^ 




(A. 79) 


(A.80) 


We choose a k such that Ck is not empty. If k < r(^o) then all A^ are equal to zero 
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therefore: 


min = min > M^± . (A.81) 

£eCfe ^ teCfc 

Let us consider k > t(^o)- Since the loop has to close, Ck contains an even number of lines. 
We separate two possibilities again: the number of lines in with (—1)'® = 1 is either 
even or odd. 


If the number of lines in Ck with (—1)^ = 1 is even, since all states propagating between 
^0 and (^1 must have (—1)'® = —1, there must be a line Iq in Lk{'H) — C with (—1)^ = —1. 
Since the proton is the lightest particle with (—1)^ = —1 then surely > Mp. We also 
pick two lines ii and £2 in Ck and we observe that ^> Mj^±. We assign 

A4 = min{£;,Mp} , 

= A^j = max 

A^ = 0, for ££ Lk{n)-{io,C,i2} . (A.82) 


E-M„ 


By using the fact that E < one can easily check that the constraints in eq. (A.66) are 
satisfied. By replacing all masses in the strange loop with M^± we get 


min \ Mf 

l^Ck ^ 


'Mk± 






if S < Mp 
if S > Mp 


(A.83) 


If the number of lines in Ck with (—1)^ = 1 is odd, then there is at least a line £i in Ck 
with (—1)'® = —1. Since the A*^ is the lightest strange particle with (—1)^ = —1 then 
surely M^^ > M^o. We also pick another line £2 in Ck- We assign 


A^^ = min 

A £2 = max 
A£ = 0 , 


\ + Ml, - Ml, 

1®’-2E- 

f„ - Ml + Ml, 

\ ’ 2E 

ioT l^Lk{n)-{£i,£2} ■ 



(A.84) 


By using the fact that E < M^- one can easily check that the constraints in eq. (A.66) are 
satisfied. By replacing the mass M^j with M^o and all other masses in the strange loop 
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with Mj^± we get 


min \ Mj — A| 

^eCfe ^ 


> 


'Mk± 



Ml,- 



V ) 

\ 


\ / 


1/2 


if ii; < 

if ii; > ^Ml, - M|.± ■ 

(A.85) 


By combining all discussed cases, we get that it is always possible to choose the quantities 
in such a way that the constraints in eq. (A.66) are satisfied and 


M{E,A) >M{E) , 


(A.86) 


where Ai{E) is defined in eq.(A.73). This concludes the proof of the theorem. 


□ 


Large-volume behaviour of the spectral density. One can easily reproduce the whole 
construction presented in this appendix around the external vertex instead of .^o- The 
general structure of the spectral density is: 




n 

eeLi 


d^pi 


p£)x 


{27r)^2Ee 








(A.87) 


where K^ g .^{E, p,ii) and n) are Euclidean kernels that include a flavour- 

violating strange loop connected to nnd respectively. The Euclidean kernels satisfy 


ln|ifJg^^(E;,p,n)| < -LM{E) + 0{lnL) , (A.88) 

thanks to theorems A.4 and A.7. The desired eq. (A.3) follows by observing that Rg ^{E, p, n) 
has a finite large-L limit and the phase-space integral generates at most powers in the vol¬ 
ume. 


B Corrections to charged-hadron masses in finite volume 


In this appendix we want to calculate the power-like finite-volume corrections to the masses 
of stable hadrons due to electromagnetic interactions at order e^. 
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Stable hadrons are identified by their (finite-volume) flavour numbers. For instance the 
charged pion has baryon number B = 0 mod 2, strangeness Fg = 0 mod 2 and electric 
charge Q = 1 mod 2. Given some flavour sector defining the target stable hadron h and 
some momentum p, we denote by \h{p),a) the lightest eigenstates of the QCD Hamiltonian 
Hq in the given flavour sector, with momentum p and with energy Ehfl{p, L). We assume 
that states with zero momentum are lighter than the others, and we refer to their energy 
as the 0{e^) finite-volume mass 

Ehfiip,L) > Eh,o{0,L) = mo{L) , if p / 0 . (B.l) 

The states \h{0),a) with a = 1,... ,ds transform under some (possibly spinorial) represen¬ 
tation of the cubic group. Normalization is chosen such that one recovers the relativistic 
normalization in infinite volume 

{h{p),a\h{p'),a') = 2Eh,o{p, . (B.2) 


The mass shift due to electromagnetic interactions is given by the Cottingham formula [27], 
which can be generalized easily to the case of finite volume by replacing the momentum 
integrals with the appropriate sums, yielding for the mass of the hadron h in finite volume 


m{L) = mo{L) 


1 ^ /■ dkoTfj,p{t,L) 

4mo{L)L\^ J 27r 


Tp,{k;L) = j d^xe-^^^h{0)\T{jp{x)UO)}\hiO)),, 


where the subscript c stands for 


(B.3) 

(B.4) 




(B.5) 


and jp{x) is the Heisenberg electric-current operator in Euclidean spacetime 

, (B.6) 

evolved with the QCD Hamiltonian Hq, and normalized in such a way that the electric 
charge is 

Q{xq) = -^ j d^xjo{xo,x) . (B.7) 

Notice that the electric current jp{x) is C-odd and therefore it is also anti-periodic in space. 
It follows that the spatial momentum k in eq. (B.4) must belong to the set n_. In eq. (B.4) 
we have also used the shorthand notation 

{h{0)\A\h{0)) = ^ ^(/i(0), a\A\h{0),a) (B.8) 

(7 


- 42 - 




which is particularly useful as the spin will play no special role in the calculation of this 
appendix. 


Because of the exponentially-raising operator in eq. (B.6), it is not obvious that the 
XQ-integral in eq. (B.4) converges. However the states that propagate in between the two 
currents have the same flavour numbers as the external state and they are odd under charge 
conjugation. In particular it follows that they cannot have zero momentum and they are 
therefore strictly heavier than the external state. In this situation the integral is shown to 
converge and an explicit calculation yields 


Tfj_u{k', L) —M^u(k] L) -\- L) , 

M^y{k-L)= j d'^xe{xo)e~"^^{h{0)\j^{x)j^{0)\h{0)). 


={h{0)\j^{0) 


Ho - mo{L) + iko 


Mo)\h{o)) 


- 2mo{L)L^ j.(0)|0) 

Hq + iko 


(B.9) 


(B.IO) 


The retarded function M^y{k', L) is analytical for any complex ko except for the simple poles 
along the positive imaginary axis. In order to make contact with the original Minkowskian 
Cottingham formula [27], the reader can easily check that T^i,{iko — e, k; oo) is the forward 
Compton amplitude for the scattering of a virtual photon with quadrimomentum k from 
the hadron h at rest. 


Formula (B.3) contains UV divergences which need to be subtracted. Notice that the electric 
current j^{x) does not require renormalization. Therefore the purely-QCD expectation 
value {h{0)\T{j^{x)ju{0)}\h{0))c appearing in eq. (B.4) is UV-hnite for x ^ 0. The operator 
product expansion of j^{x)ju{0) implies that T^u{k, L) vanishes like k~'^ (up to logarithms) 
at large k, which makes the integral in eq. (B.3) logarithmically divergent. Following [28] 
we renormalize the Euclidean Cottingham formula by introducing a Pauli-Villard regulator 
for the photon propagator and by adding appropriate counterterms: 


m{L) = mo{L)+ 


-|- lim < — 

A^OO 


4mo(T) L3 


E 

ken_ 


dkoT^fj,{k;L)A‘^ 

2tt k‘^{k? + h?) 


+ (h(0)|C(A)|h(0))A . (B.ll) 


In this formula we assume that the regulator needed to dehne T^i,[k,L) has been already 
removed. The counterterms have the form: 


C'(A) = C6i(A)( 9^^(0) -7 ^c/(A)m/V’/V’/(0) , 

f 


(B.12) 


where is the (Euclidean) energy-momentum tensor. Since QCD-|-QEDc is a local the¬ 
ory, the coefficients c(A) can be chosen to be L-independent by choosing renormalization 
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conditions in infinite volume. 


We are ready now to manipulate eq. (B.ll) in order to extract the power-like finite-volume 
corrections to the mass. 

Lemma B.l. The QCD quantities appearing in eq. (B.ll) have only exponentially-suppressed 
finite-volume correetions, 

mo{L) - mo(oo) = , 

L) - oo) = 0{e~^^^^) , for any real A: / 0 , 

(h(0)|C(A)|h(0))c - lim (h(0)|C(A)|h(0))c = . (B.13) 

L^oo 


Proof. A possible proof of this lemma, which we will not give here, can be obtained by using 
intermediate results and theorems in [22], under the assumption that the leading finite- 
volume corrections are described by some arbitrarily-complicated Lagrangian massive field 
theory with small couplings, which effectively describe the dynamics of hadrons at large 
distance. All above quantities can be decomposed in terms of dressed propagators and (IPI) 
proper vertices with possible insertions and with two on-shell external legs. For all these 
quantities, the general conclusions of theorems 2.4, 2.5 and 2.6 hold, leading to a proof of 
the lemma. Notice that for a general theory, the finite-volume effects on the masses are 
however this is not the case in QCD [29]. Some of the technology of [22] is 
adapted to the case of C*-boundary conditions in app. A. □ 


Thanks to lemma B.l, we can write for the finite-volume correction to the mass 


Am{L) = m{L) — m(oo) 


1 


lim < —r ^ 


4mo A->-oo I Lfi 


d^k [ fdko T^^{k)A^ 
( 27 r )3 I y 27 r k‘^{k‘^ -\- 


+ 


ken_ 

where it is understood that we mean L = oo whenever we drop the L dependence. 


We introduce an arbitrary function r]{z) of a real variable z with the properties: (a) r]{z) is 
infinitely differentiable for any value of z, (b) rj{z) = r]{—z), (e) r]{z) = 1 for \z\ < M^/2 for 
some arbitrary M > 0, (d) r]{z) = 0 for \z\ > M^. We rewrite the finite-volume correction 
to the mass as 


Am{L) 


4mo 


1 

U 


E 

ken_ 


d^k 

(27r)^ 


r](k^) 


dkp Tfj,^{k) 
2tt 


+ R{L) 


(B.15) 
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where the reminder R{L) is 


R{L) = 


lim 

A^OO 


ken_ 


e 

4mo 


[l-r/(k2)] j 


^}TA(k) + 0(.- 

'dko T^^{k)A^ 

27r k‘^{k‘^ + A^) 


) , 


(B.16) 


Lemma B.2. The infinite-volume Euclidean amplitude Tfj^fj,{k) is infinitely differentiable 
for any k G M^/{0}. 


Proof. See app. B.l. 


□ 


Thanks to lemma B.2 and to the fact that the factor 1 — ?/(k^) regularizes the singularity 
in k = 0, XA(k^) is inhnitely differentiable in k G The reminder R{L) is the difference 
between the integral of a smooth function and its approximation as a Riemann sum, which 
vanishes in the inhnite-volume limit faster than any inverse power of L, 

lim L‘^R{L) = 0 , for all a; > 0 . (B.17) 

Tj—¥ CXD 


By plugging the eq. (B.9) into eq. (B.15), and by folding the feo integral we get 


Am{L) 


e 


2 


2mo 


1 

L3 


E- 

ken_ 


d^k ^ 
(27r)3 


ri(k^) 


dkp Mf^ffk) 
271 k"^ 


+ R{L) . 


(B.18) 


The integrand is holomorphic in the half plane Im/co < 0 except for the single pole in 
ko = —^|k|. Therefore the ko integral can be calculated as a Cauchy integral by closing the 
contour at inhnity in the lower half plane yielding 


4m„ i L3 E_ / (2^)3 




(B.19) 


We will see that the power-law hnite-volume corrections come from the behaviour of the 
integrand in eq. (B.19) around k = 0. 


Lemma B.3. Because of rotational symmetry, the on-shell retarded function M^^(—z|k|,k) 
is a function o/k only via |k|. R can be decomposed as 


M^^(-z|k|,k)||^l^^ 


M-i 

K 


-P M{k) , 


(B.20) 
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where the function A4 (k) is analytical for complex values of k in a neighbourhood of zero, 
and 


M-i = -2moq^ , 

where q is the electric charge of the hadron h. 


(B.21) 


Proof. See app. B.l. 


□ 


We can use now Poisson summation formula in order to express the discrete sum in 
eq. (B.19) over spatial momenta in terms of Fourier integrals. Since the momentum k 
belongs to the anti-periodic set n_, an extra sign appears in Poisson summation formula 

E /O') = E (-!)<“’/ /(k) , (B.22) 

ken„ neza '' '' 


where (n) has been defined in eq. (2.22). The term n = 0 in the previous expression cor¬ 
responds to the infinite-volume integral. By plugging the definition (B.20) into eq. (B.19), 
and after calculating the angular integral in k, we get 


Am{L) = — 


e^M-i ^ (—1)^"^ j 2\Sin(K|n|L) 


SmoTT^L 




dK r]{K'^ 


neZa/lo} 


Svr^moT 


—I—I— J dn sm{K\'n\L) P R{L) . 

(B.23) 


neZa/{0} 


We exploit the arbitrariness we have in choosing the function ri{K^) and assume that it has 
support in the analyticity domain of Thanks to lemma B.3 the function rj{K?‘)Ai[iP) 

is smooth for any k > 0, and has all right derivatives in k = 0. The expansion in powers of 
1/L can be written in terms of the following generalized zeta function 

«>)= E (b.24) 

neZa/lo} 

which is analytically extended to a meromorphic function in the whole complex plane, and 
holomorphic for Re s > 0. 


The first integral in (B.23) can be understood by defining the function 


V{x) = J ^ 


(B.25) 


Since r](k‘^) is a Schwartz function so is ii(x), and in particular it decays at infinity faster 
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than any inverse power of x. The sum 


2 (-!)<") 

TT Ini 

neZs/jO} 



dn r]{K^) 


sin(K|n|L) 

K 


—I—I— / dx f}{x) (B.26) 
neZs/fO} J-\n\L 


converges to ^(1), and the corrections decay faster than any power in 1/L. The second inte¬ 
gral in (B.23) has a Taylor expansion in (|n|L)“^ that can be extracted by using iteratively 
the identity 


f 


dn /{k) sin(xK) 


1 1 

^/(O)-^ / dK f”{K)sm{xK) . 

X x^ Jq 


(B.27) 


Putting everything together we obtain the desired expansion of the finite-volume corrections 
to the mass 


Am{L) 


16mo7rL 


e(l) 


Svr^mo 


OO 


E 


(-1)' 

£ 2 + 2 £ 


+ 2 £) . , 


(B.28) 


where the dots stand for contributions that decay faster than any power of 1/L, and A42£ 
is the (2^)-th derivative of JA{k) in k = 0. Notice that the on-shell forward Compton 
amplitude is given by 


r(k2) = T^^(z|k| - e,k) = M^^(z|k|,k) +M^^H|k|,k) =Af(|k|) + Af(-|k|) , 


(B.29) 


therefore the coefficients ^-re trivially related to the derivatives 

(2f-)! d^ {2£)l 

2(£!) d(K2y^^^^ ~ 2{£\)^^ 


(B.30) 


of the on-shell forward Compton amplitude for the scattering of soft photons on the hadron 
h at rest. 


We also notice that the coefficients Af_i and To depend only on the mass and charge of 
the hadron, and not on its spin or internal structure. For the scattering amplitude we use 
the classical result [30, 31] (also reviewed in section 13.5 of [32]): 

To = lim ^^;,(^|k| - e, k) = -Aq^ . (B.31) 


We conclude this appendix by providing a representation of the zeta function ,^(s) defined 
in eq. (B.24) which is useful for numerical calculation. We use the identity 


1 


1 


n 


*/2 r(s/2) 




1 


r(s/ 2 ) 




(B.32) 
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we plug it into eq. (B.24) and we use the Poisson summation formula 



3 


^-wn?+i-K{n) 

neZs 


ken_ 


Au 


(B.33) 


only in the integral over u G [0, u*]. At this point all integrals can be calculated explicitly 
in terms of the upper incomplete gamma functions: 




The upper incomplete gamma function r(T, z) is defined for all complex values of r except 


non-positive integers, and it decays exponentially as \z\ —)• oo. Therefore the infinite sums in 
the previous formula are rapidly convergent. Also this representation is valid for all values 


of s needed in the mass formula. The splitting variable rt* > 0 is completely arbitrary and 


can be used to check the result of numerical calculation. 

B.l Analyticity properties 

In this subsection we work in Minkowski spacetime with metric g = diag{l, —1, —1, —1). We 
also set L = oo. We introduce the Minkowskian electric current: 

J^{x) = ed^o/^o-xP) j^(Q)g-*(xoHo-xP) ^ 35^ 

which is related to the Euclidean one introduced in eq. (B.4) via 


J\0) =-ijoiO) , jH0)=jki0). 


(B.36) 


While in finite volume, because of C*-boundary conditions, eigenstates of the momentum 
are also eigenstates of the charge-conjugation operator, this is not necessarily true in infinite 
volume. We perform a change of basis which does not affect the quantities we are interested 
in, and we choose to work with simultaneous eigenstates of energy, momentum and electric 
charge. 

We consider the retarded two-point function in the forward limit 



(B.37) 


which is related to the function introduced in eq. (B.IO) in infinite volume through 
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a Wick rotation 


M^^(/co,k) = -W+(-^feo,k) . (B.38) 

The subtraction of the disconnected part in eq. (B.37) can be expanded: 

I d4xe*'="0(xo)(h(p)|J^(x)J^(O)|h(O))e 

= I d^xe*'^"0(xo)|(h(p)|4(x)J^(O)|h(O))-2i7(p)(27r)3h3(p)(o|4(x)J'^(O)|O)| . 

(B.39) 

Notice that in this formula we cannot just take p = 0 because the disconnected contribution 
gives a geometrical divergence proportional to h^(0). Therefore the limit in eq. (B.37) is 
essential in order to define properly the subtraction. However notice that for any p 7 ^^ 0 the 
delta function vanishes exactly and it does not contribute to the limit (limp_>o = 0 ), 
allowing us to write equivalently 

W+{k) = i lim f e*^^0(xo)(/i(p)| Jp(x) J^(0)|/i(0)) , (B.40) 

p-^oj 

It is possible to prove that this limit is finite, which we will assume in the remaining of this 
appendix. ‘ 

By calculating the coordinate integral in eq. (B.40) one gets 

W+{k) = lim (ft(p)l J,(0) ~ P ~ J'‘(0)|ft(0)) , (B.41) 

p-s>0 Hq — h[p) — Kq — l€ 

where we have introduced the energy of the external states 

E{p) = + p2 . (B.42) 

Notice that, since the external state corresponds to a stable hadron, the states propagating 
in between the two currents in eq. (B.41) have energy not smaller than mo. Therefore the 
retarded function W+(fe) has poles only for non-negative values of /cq. 

It is useful to separate the single particle component from the continuous part of the spec- 


^ By means of the LSZ reduction formula, the function C+{k) dehned in (B.50) can be expressed as 
linear combinations of the reduced Green’s functions dehned in eq. (16.52), chapter 16 of [33]. On the other 
hand, as we will notice later on, W+{k) is uniquely determined by C+{k). The hniteness of the limit p —>■ 0 
in eq. (B.50) and consequently in eq. (B.40) derives from the analyticity properties of the reduced Green’s 
functions stated in Theorem 16.8, chapter 16 of [33]. The reader should notice that the connected part and 
the limit are systematically dropped in the classical literature, e.g. [27]. 
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trum: 


W+{k) = + ^Mp(fco,k^) , (B.43) 

£/(k) — mo — kq-\- le 

1 

= 4 E ^n^^\ih{0),cT\MO)\h{k),a’)\^ . (B.44) 

fX^C^Cj' 


We want to study the analyticity properties of Zip(k^) and Zmp(A:o, k^) in the spatial 
momentum k (the analyticity properties in /cq are obvious from the spectral decomposi¬ 
tion (B.41)), which we summarize here: 

1 . Zip(k^) is analytical for any real value of k, and can be analytically continued to a 
complex neighbourhood of k^ = 0 ; 

2. The Euclidean function Zmp(—*A: o, k^) is analytical for any real value of k; 

3. The on-shell function Zmp(—| k|, k^) is analytical for any real value of k; 

4. The on-shell function ZMp(|k|, k^) is analytical for real values of k in a neighbourhood 
of k = 0, and can be analytically continued to a complex neighbourhood of |k| = 0 . 

From these properties it follows that: 

1. The off-shell Euclidean Compton amplitude 


r(/co, k) = M^^{ko, k) -h Mfj_^{-ko, -k) 


2[E{k) - mo]Zip(k^ 
[E;(k) - mo ]2 -h fcg 


- Zmp(*A:o, k^) — Zi/ipi—iko, k^) 


(B.45) 


is analytical for any real value of k 7 ^ 0 (lemma B.2). 

2. The on-shell quantity 


M^^(-i|k|,k) = - 


Zip(k2) 


E{k.) — niQ + |k| 


— ■^Mp(—|k|,k^) , 


(B.46) 


as a function of |k|, admits a meromorphic extension to a complex neighbourhood of 
|k| =0. In particular it admits a Laurent series in |k| = 0, the first term being: 


M^^(-i|k|,k) = -%M + o(|k|0) 


(B.47) 


We will show that, eq. (B.65), Zip(O) = 2moq^ (lemma B.3). 

3. The on-shell Compton amplitude 


T(z|k| — e, k) = — 


2 [E;(k) - mo]Zip(k^) 
[-E(k) — mo]^ — k^ — ze 


— ZMp(|k|, k^) — Zmp(—| k|, k^ 


(B.48) 


is an analytic function of k^ in a complex neighbourhood of k^ = 0. This follows from 
the analyticity properties discussed above, from the fact that the odd powers in |k| 
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generated by the expansion of ZMp(±|k|, k^) cancel ont, and from the fact that 


2[E{k) - mo] 

k-5>o [-E(k) — mo]^ — k^ — ze 


1 

niQ 


(B.49) 


The full analyticity properties of and consequently of Zip(k^) and Zyipiko, k^), can 

be derived by the analyticity properties of four-point reduced Green’s functions discussed 
in chapter 16 of [33]. However we provide here a hopefully more digestible proof of the 
particular properties we are interested in, based on the Jost-Lehmann-Dyson representation 
of the expectation values of certain retarded commutators. We also point out that the same 
analyticity properties we are interested in can also be obtained by assuming an effective 
theory describing hadrons and by using results and methods discussed in section 2.4 of 
ref. [22] and in appendix A. 


Analysis of Zip(k^). We notice first that Zip(k^) can be extracted also from the retarded 
commutator 

C+{k) =i lim f d^xe*^’^6'(xo)(/i(p)| [J^(x), J^(0)] |/i(0)) 

=W+{k)+ W-{k) , (B.50) 

where W-{k) is a functions with poles only for negative value of /cq. Therefore the following 


reduction formula holds 

lim [^'(k) — niQ — ko]C+{k) 

ko^E{U.)—mo 

= lim [iil(k) — mo — fco]lT+(A;) = Zip(k^) . (B.51) 

ko^E(k.)—mo 

Then we extract the single-hadron pole from both orderings of the retarded commutator 
by means of the following trick. We introduce the auxiliary retarded commutator 

C'+(fc) = z lim [ d^x e^^^9{xo){h{p)\[Jfj_{x), J^{0)]\h{0)) , (B.52) 

Jfj.{x) = (—□-|-2zmo9o) J^(x) . (B.53) 

The relation between this retarded commutator and the original one is obtained through 
integration by parts of the differential operator (—□ -|- 2zmo5o) 

(/c^ -|- 2moko){k'^ — 2moko)C+{k) = C+{k) + P{k) . (B.54) 

The boundary term has the form 

P{k) = i lim f Px6{xo)e^^^ {h{p)\[T>Jn{x),J'^{0)]\h(0)) (B.55) 

p^Oj 


where P is some local differential operator. The integrand of P{k) involves only commuta¬ 
tors of local operators at equal time, which are linear combinations of delta functions and 
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their derivatives. Therefore P{k) is a polynomial in the quadrimomentum k. In terms of 
the auxiliary retarded commutator, the reduction formula reads 


Zip(k 2 ) 


1 

8 moi?(k)[iil(k) — mo] 


lim [C+ik)+P{k)]. 

ko^E[k.)—mo 


(B.56) 


The analyticity properties of the modified retarded commutator can be exposed by means 
of the Jost-Lehmann-Dyson (JLD) representation [34, 35], 


C+{k) 



(Pud\^ p{u, 
up - + i{ko - uo)e 


(B.57) 


where the JLD spectral function p{u, A^) is uniquely determined by the retarded commu¬ 
tator. The integration domain S encodes all known information about the spectrum, and 
can be represented as the set of 5-tuples {u,}?) such that 


V^(u - k)2 -I- A^ + uo> JMl -I- k^ — m 

< ,_ \ - , (B.58) 

\/(u — k)2 -|- A^ — uq> M| -|- k^ — m 

for any value of the momentum k. The masses Mi and M 2 are determined in the following 
way. Consider the commutator 


lim [pxP^ph{p)\[j^{x),j>^m\m) 

= ^im(^(p)l«4‘(0)(2vr)^[<5(P - p - k) - (5(T’ -p + fe)]p=(£;(p),p) J^(0)]|/i(0)) , (B.59) 

with P = The two delta functions come from the two different orderings of the 

currents in the commutator. Mi and M 2 are the masses of the lightest states propagating in 
between the two currents in the first and second ordering respectively. Had we considered 
the original current Jp(x), the lightest state would have been the hadron h itself. However 
it is easy to check that the insertion of the operators (—□ -|- 2^mo5o) kills the contribution 
of the single-hadron states in the above commutator, therefore 


Ml = M 2 = m -|- A 

where A > 0 is some mass gap (if no bound states exist A = 2mp. 


(B.60) 


The relevant limit for the reduction formula (B.56) is 


/(k2)= hm [C+{k)+Pik)] 

ko^E{K)—mQ 

=P{E{k),k) + 


(PudX^ p{u, A^ 


[£'(k) — mo — uoP — (k — u)^ — A^ 


(B.61) 
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The denominator vanishes only if 


uq ± \/(k — u)2 + = -E'(k) — mo , (B.62) 

which is satished for no real value of k if [u, A^) is in the domain S. 

As pointed out in [35], if {u, A^) is in the domain S then necessarily 

l^^ol + |u| < mo . (B.63) 

Thanks to this, it is easy to show that two positive constants k and a exist such that the 
denominator in eq. (B.61) is limited from below by 

I[i?(k) — rriQ — tio]^ — (k — u)^ — A^| > |tto — — A^| — ajkl > — Q:|k| (B.64) 

for any {u, A^) G S and for any complex k such that |k| < k. In the last step we have used 
eqs. (B.58) for k = 0. From the above bound it is clear that if k is small enough, then the 
denominator never vanishes. Therefore /(k^) can be continued by analyticity to a complex 
neighbourhood of k^ = 0 . 

From eq. (B.56) it might seem that Zip(k^) has a singularity in k —0. However this limit 
is fixed by symmetries: 

1 _ 

^im Zip(k2) = ^iin — ^ (h(0), cj| J^(0)|/i(k), ct')^ = 2mog^ , (B.65) 

11 ,( 7 ,< 7 ^ 

where q is the electric charge of the hadron h. This relation implies that /(O) = 0 and 
Zip(k^) is an analytic function for any real value of k^, and for complex values of k^ in a 
neighbourhood of zero. 


Analysis of Zmp(^ 0;k^). ZMp(A;o,k^) is obtained by selecting all poles in IF+(fc), or 
equivalently in C+{h), with Refco > E(k) — mo. In this case we find more convenient to 
write C^{k) in terms of the auxiliary retarded commutator 

C+{k) = t lirn^ f d'^x e^^^6{xo){h{p)\[J^{x), J^{0)]\h{0)) . (B.66) 

p->o ^ j 

In complete analogy to eq. (B.54), the original retarded commutator can be written in terms 
of the auxiliary one as 

(A:2 + 2moA:o)C+(A:) = C+{k) + P{k) , (B.67) 

where P{k) is a polynomial in the quadrimomentum k. The above equation can be inverted 
by noticing that all poles of C+(x) have negative imaginary part. We introduce a JLD 
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representation for the retarded commutator Cj^{k) and we get 


d'^udX^ p{u, A^) 


1^2 _j_ 2 moko + i{ko + mo)e _ ^2 + ^(A;o - uo)e 


(B.68) 


The integration domain S is the set of 5-tuples (u, A^) such that 

J -^/(u — k)2 -|- A2 -|- Mo > ^/{m + A)2 -|- — m 

\ -^/(u — k)2 -|- A2 — uq > \/m^ + k'^ — m 


(B.69) 


for any value of the momentum k. The denominator outside of the integral in eq. (B.68) 
has poles for Re/co < — rriQ which do not contribute to ZMp(A:o,k^). The integrand 

can be decomposed in partial fractions 


1 


{k — m)2 — A2 -|- i[kQ — Mo)e 

X = ^/{k- U)2 + A2 , 


— f - - - 

2X \ko — uo — X + te 


1 


ko - Uq + X + ie 


(B.70) 

(B.71) 


and only the first one contributes with a pole to ZMp(A:o,k^). By calculating the residue 
at this pole we get 


■^Mp(^ 0 )k^)= / d^udX‘^ 


p{u, A^ 


2X[(mo + a: + mo)2 - E{ky] ko-uo-X + ie 


(B.72) 


Using the definition of the domain S it is straightforward to check that the denominator 
2 X[(mo + X + m-o)^ — ill(k)^] never vanishes for any real value of k and for any value of 
(m, A^) G 5. We derive some particular properties. 


The Wick-rotated function 


ZMp{—iko,k^) = / d^udX^ 


p{u,X^) 


(B.73) 


2A:[(mo + a: + mo)2 - E{ky] -zko -uo-X 
is analytical for any real value of k, since the two denominators never vanish (as uq+X > 0). 


The on-shell function 


■^Mp(—|k|, k^) 



_ p{u,X^) 

2X[{uQ + X + moY 


1 

.E;(k)2] -|k| - Mo - a: ■ 


(B.74) 


is analytical for any real value of k, since the two denominators never vanish (as |k| -|- 
Mo -|- W > Mo -|- Ai >0). Moreover, as a function of |k|, Zmp(—| k|,k^) can be analytically 
continued to a complex neighbourhood of |k| = 0. 
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The on-shell function 



P{u,\^) 

2X[{uo + X + mof-E{\if] 


1 

k| — uo — ^ + ^e 


(B.75) 


is analytical for real values of k such that 


|k| < A , (B.76) 

as in this range the denominator |k| — uq — X can be shown not to vanish for any value of 
{u,}?) (z S. As a function of |k|, ZMp(|k|,k^) can be analytically continued to a complex 
neighbourhood of |k| = 0. 


C Classical vacua of compact QEDc 


We consider an abelian gauge field on a lattice with C* boundary conditions along the 
directions included in the set C 


U{x + L^,p) 


U{x, p) if p ^ C 
U{x, p)* if /r G C 


(C.l) 


where the coordinates are integer numbers in the range 


We assume direction p = 3 C*-periodic, and direction p = 0 periodic. 


(C.2) 


We want to characterise all gauge-field conhgurations corresponding to absolute minima of 
the Wilson action. In terms of the plaquette P{x, p, v) the minimum condition reads 

P{x, Pjv) = 1 . (C.3) 

We can always gauge-transform to axial gauge along a given direction p, i.e. to a gauge in 
which all the link variables U{x,p) are equal to one except the ones on the hyperplane vr^ 
defined by the equation 

TTfj, : x^ = Lf, - 1 . (C.4) 

Because of condition (C.3) it is easy to show that we can gauge-transform to simultaneous 
axial gauge for all directions. We will refer to those link variables that are different from 
unity as active link variables (see figure 7). 

Plaquettes at the intersection of two distinct vr-planes involve four active link variables, all 
other plaquettes on the vr-planes involve two parallel link variables. Constraint (C.3) on the 
latter ones, together with the fact that unity is left unchanged by the boundary conditions. 
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Figure 7: A two-dimensional representation of the problem discussed in this appendix. 
Direction = 0 is periodic, while direction /r = 3 is C*-periodic. In simultaneous axial 
gauge, the only links that are different from unity are the ones represented with a thick line 
and an arrow {active link variables). The condition that the blue plaquette be equal to one 
implies that the two active link variables in the plaquette are equal. Nontrivial constraints 
come from minimum condition for the red plaquette at the intersection of the vro and tts 
hyperplanes. 


implies 

U{x, ^) = U{x + 0, n) , for any x G vr^, u ^ fa . (C.5) 

Using this equation recursively we get that, given some direction fi, all active link variables 
along fi are equal to each other. We will define 

= U{x,fi) , for any x G vr^ . (C.6) 

We use now the minimum condition (C.3) for plaquettes at the intersection of two distinct 
TT-planes. Let us consider first the plaquette in some point x G TT^nvra where ^ is a periodic 
direction (see figure 7): 

1 = P(x, fi, 3) = , (C.7) 

which implies 

Wf, = ±i, afi^c . (C.8) 
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If ;U 7 ^ 3 is a C* direction we get instead 


1 = P{x,^i,2,) = , 


which implies 


= ±11^3 , if M G C • 


(C.IO) 


Finally we show that IF 3 can be set to 1 with a gauge transformation. Let re be a complex 
number such that IF 3 = w~‘^ and we define the gauge transformation 


A(x) = rc , for 0 < — 1 , 


(C.ll) 


and extended outside the above domain by means of the boundary conditions 


A(x + Lfj) — 


A(x) \i jj, ^ C 

A(x)* if /U G C 


(C.12) 


First notice that this gauge transformation preserves the gauge-held boundary conditions 
and the axial gauge. All active link variables along periodic directions are left unchanged 
under this gauge transformation. If ^ is a C*-direction, the active link variable along ^ 
transforms like 


Wf, ^ wW^w = LFg-^LF^ = ±1 , 

with the particular case of 

IF3 ^ wW^iW = IFg-^IFa = 1 . 

This concludes the proof of part 1 of the following proposition. 


(C.13) 


(C.14) 


Proposition C.l. Let U{x,fi) be a gauge eonfiguration that minimizes the Wilson aetion. 
1. A veetor z satisfying the eonditions 


2:3 = 1 , = 1 


(C.15) 


exists sueh that U{x,fi) is gauge-equivalent to the gauge eonfiguration 


Ufix,^) = 


Zfi if x^ = L^-l 
1 otherwise 


(C.16) 


2. The veetor z is unique. 


Uniqueness is proven by noticing that the vector 2 : is therefore uniquely determined by the 
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original gauge configuration U{x,n) 


fVF(M) ifM0C 

\w{^)W{ 3 )-^ if^uGC 


where we have introduced the Wilson lines 




s=0 


(C.17) 


(C.18) 


It is easy to show that W{fi) is gauge invariant if and only if /x is a periodic direction, 
while the L-shaped parallel transport W{^)W{3)~^ is gauge invariant if and only if /x is a 
C*-direction. 


D Anatomy of the sign problem 


Integration of the fermion fields in a periodic setup yields the determinant of the Dirac 
operator. This result relies on the fact that the Grassman variables ^{x) and x/)(x) are 
independent, which is not true in the case of C* boundary conditions. A possible way to 
get an explicit expression for the fermionic path integral is to use the change of variable 


'ip±{x) 


V’(x) ± C 

71 


(D.l) 


and to define the new two-component field 


r]{x) = 


y-l^P-{x) j 


(D.2) 


It is straightforward to verify that C* boundary conditions for the field V'(x) are equivalent 
to 


r]{x + Lk) = Kr]{x) , 



(D.3) 


where the matrix K acts on the two components of r/. We will refer to these boundary 
conditions as K boundary conditions. 


By using the identity for the Wilson-Dirac operator (valid for a general non-abelian gauge 
theory) 

C-^D[VfC = D[V*] , (D.4) 
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and a few lines of algebra, one can write the fermionic action in terms of the new fields 


SF = i^D[V]2l; = -^rj^CD[J{V)]T] , (D.5) 

where Dj = D[J'{y)] is the Wilson-Dirac operator calculated with the gauge field Jiy) 
defined as 


J{V) = I 2 < 8 ) Re y + J (g) Im y , 



(D. 6 ) 


The matrices I 2 and J act on the two components of rj. Notice that JiV) defines a 
representation of the gauge group, unitarily equivalent to the representation defined by V. 
Integration of the fermionic action in the form obtained in eq. (D.5) yields 

[ = [ Vrj = Pf^ CDj , (D.7) 

JC* b.c.s Jk b.c.s 

where the subscript K reminds that the derivative appearing in the Dirac operator are de¬ 
fined on the space of fields satisfying K boundary conditions, and CDj is an antisymmetric 
complex matrix. In eq. (D.7) Pfi^- CDj is the Pfafhan of CDj that, by using the algebraic 
identities 


(Pfx CDjf = Det^ CDj = Det^ Dj , (D. 8 ) 

can be related to the determinant oi Dj. Algorithms for the lattice simulation of theo¬ 
ries involving Pfafhans have been discussed in the context of C* boundary conditions or 
the closely-related G-parity boundary conditions and also in the context of lattice super- 
symmetric models (see [36-39] for a list of references on this subject). 


We shall now discuss if a sign problem is associated to Pfi^ CDj. By using eq. (D. 8 ) and 
the 75 -hermiticity of the Dirac operator, one concludes easily that the squared Pfafhan is 
real. We want to show now that a stronger result holds: the Pfafhan itself is real. Let us 
consider the Pfafhan of the auxiliary operator C{Dj — s) for a generic complex number s. 
This Pfafhan is a polynomial in the matrix elements and in particular in s, 

Pfi^ C{Dj -s) = n(s - A„)™“ , (D.9) 

a 

where the Aq’s are distinct roots. The overall normalization is determined by the value 
of the Pfafhan in the s —)• 00 limit. By using the relation between the Pfafhan and the 
determinant we calculate the characteristic polynomial oi Dj 

BetK {Dj -s) = [PfK C{Dj - s)]2 = His - A„)2™“ . (D.IO) 

a 

The Aa’s are the roots of the characteristic polynomial of Dj, i.e. they are the eigenvalues 
oi Dj. Notice that the algebraic multiplicity of A^ is 2ma- Because of 75 -hermiticity either 
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the eigenvalues of Dj are real or they appear in pairs of complex conjugates. Since all 
multiplicities are even, the determinant is positive if s is real, and consequently the Pfafhan 
is real. For s = 0 one gets 

PfKCDj= n |A^“|2. (D.ll) 

Q:|ImAcK=0 a|ImAQ;>0 


Once established that the fermionic Pfafhan (D.7) is real, we need to wonder about its sign. 
From eq. (D.ll), clearly the Pfafhan is negative only if the Dirac operator Dj has some 
negative eigenvalues, which can happen with Wilson fermions. However, in the continuum 
limit, the real part of the eigenvalues of the Dirac operator is always positive (and equal to 
m) therefore the Pfafhan is positive. At hnite lattice spacing the fermionic Pfafhan (D.7) 
has a mild sign problem that is completely analogous to the single-havour case with periodic 
boundary conditions. 
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